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SPECIAL NOTES 


1. API PUBLICATIONS NECESSARILY ADDRESS PROBLEMS OF A GENERAL 
NATURE. WITH RESPECT TO PARTICULAR CIRCUMSTANCES, LOCAL, STATE, 
AND FEDERAL LAWS AND REGULATIONS SHOULD BE REVIEWED. 


2. APIIS NOT UNDERTAKING TO MEET THE DUTIES OF EMPLOYERS, MANU- 
FACTURERS, OR SUPPLIERS TO WARN AND PROPERLY TRAIN AND EQUIP 
THEIR EMPLOYEES, AND OTHERS EXPOSED, CONCERNING HEALTH AND 
SAFETY RISKS AND PRECAUTIONS, NOR UNDERTAKING THEIR OBLIGATIONS 
UNDER LOCAL, STATE, OR FEDERAL LAWS. 


3. INFORMATION CONCERNING SAFETY AND HEALTH RISKS AND PROPER 
PRECAUTIONS WITH RESPECT TO PARTICULAR MATERIALS AND CONDI- 
TIONS SHOULD BE OBTAINED FROM THE EMPLOYER, THE MANUFACTURER 
OR SUPPLIER OF THAT MATERIAL, OR THE MATERIAL SAFETY DATA SHEET. 


4. NOTHING CONTAINED IN ANY API PUBLICATION IS TO BE CONSTRUED AS 
GRANTING ANY RIGHT, BY IMPLICATION OR OTHERWISE, FOR THE MANU- 
FACTURE, SALE, OR USE OF ANY METHOD, APPARATUS, OR PRODUCT COV- 
ERED BY LETTERS PATENT. NEITHER SHOULD ANYTHING CONTAINED IN 
THE PUBLICATION BE CONSTRUED AS INSURING ANYONE AGAINST LIABIL- 
ITY FOR INFRINGEMENT OF LETTERS PATENT. 


5. GENERALLY, API STANDARDS ARE REVIEWED AND REVISED, REAF- 
FIRMED, OR WITHDRAWN AT LEAST EVERY FIVE YEARS. SOMETIMES A ONE- 
TIME EXTENSION OF UP TO TWO YEARS WILL BE ADDED TO THIS REVIEW 
CYCLE. THIS PUBLICATION WILL NO LONGER BE IN EFFECT FIVE YEARS AF- 
TER ITS PUBLICATION DATE AS AN OPERATIVE API STANDARD OR, WHERE 
AN EXTENSION HAS BEEN GRANTED, UPON REPUBLICATION. STATUS OF THE 
PUBLICATION CAN BE ASCERTAINED FROM THE API AUTHORING DEPART- 
MENT [TELEPHONE (202) 682-8000]. A CATALOG OF API PUBLICATIONS AND 
MATERIALS IS PUBLISHED ANNUALLY AND UPDATED QUARTERLY BY API, 
1220 L STREET, N.W., WASHINGTON, D.C. 20005. 


Copyright © 1993 American Petroleum Institute 
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FOREWORD 


API publications may be used by anyone desiring to do so. Every effort has been made 
by the Institute to assure the accuracy and reliability of the data contained in them; however, 
the Institute makes no representation, warranty, or guarantee in connection with this pub- 
lication and hereby expressly disclaims any liability or responsibility for loss or damage re- 
sulting from its use or for the violation of any federal, state, or municipal regulation with 
which this publication may conflict. 

Suggested revisions are invited and should be submitted to the director of the Measure- 
ment Coordination, Industry Services Department, American Petroleum Institute, 1220 L 
Street, N.W., Washington, D.C. 20005. 
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INTRODUCTION 


This document is the Documentation File to API Publication 2518, Second 
Edition [A7]*. 


The purpose of the Documentation File is to present detailed technical 
information related to the development of API Publication 2518 that includes: (1) 
the development of theoretical equations; (2) comparisons with test data; (3) a 
sensitivity analysis of the loss equation; and (4) other pertinent information 
that was developed during the preparation of API Publication 2518. 


The Documentation File is divided into two main parts: Sections A through H 
pertain to the standing storage loss, and Sections I through L pertain to the 
working loss. 


The standing storage loss equation in the Second Edition [A7] is different 
then that in the First Edition [A6]. Sections A through H present the 
development of the new standing storage loss equation. 


The working loss equation in the Second Edition [A7] is the same as that in 
the First Edition [A6]. Sections I through L contain development information 
that originally appeared in the First Edition. 


Section R contains a list of important References that were reviewed in 
developing the Second Edition. These references are cited in various sections of 
the Documentation File. 


is Numbers in brackets refer to the numbered references listed at the end of 


this Documentation File. 
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SECTION A 
NOMENCLATURE 

DESCRIPTION 
Vapor pressure function constant 
Vapor pressure function constant 


Vapor space expansion factor 
Number of moles 


. Pressure 


Pressure change 

Ideal gas constant (10.731) 
Temperature 

Temperature change 

Volume 

Volume change 

Mole fraction in the vapor phase 


SUBSCRIPTS 


Air 

Atmospheric 

Breather Vent 

Breather Vent Pressure Setting 
Breather Vent Vacuum Setting 
Liquid 

Liquid Average 

Tota} 

Stock Vapor 

Vapor Average 


Initial Condition or Minimum Condition 
Final Condition or Maximum Condition 
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OR 
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mole 
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psi 


psia ft3/Ibmole OR 
OR 

OR 
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ft3 

mole fraction © 
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Al1.0 INTRODUCTION 


This section of the Documentation File to API Publication 2518, Second 
Edition, contains a derivation of the equation for the vapor space expansion 
factor, Ke. This equation is derived from the ideal gas law and from the 
pressure, temperature and volume conditions that exist in the vapor space of a 
fixed-roof tank containing a volatile liquid stock. 


Section A2 presents a derivation of the vapor space volume change due to 
thermal breathing. This derivation closely follows that originally derived by 
Boardman [1]* and that, presented at the API “Symposium on Evaporation Loss’ of 
Petroleum from Storage Tanks" November 10, 1952 [2]. 


Section A3 defines the vapor space expansion factor, Ke, and develops the 
equation that may be used to calculate this factor. 


Section A4 describes various simplifications that can be made to the 
equation for the vapor space expansion factor to permit ease of calculation with 
little loss in accuracy. 


A2.0 VAPOR SPACE VOLUME CHANGE 


Figure Al is a schematic illustrating the tank vapor space thermal breathing 
process in a fixed-roof tank that is partially filled with a volatile liquid 
stock and equipped with a pressure-vacuum vent. 


During the thermal breathing process, the pressure, volume: and temperature 
vary from minimum condition 1 to maximum condition 2. At conditions 1 and 2, the 
total absolute pressure in the vapor space is P and Pz, respectively, where: 


Py = Patm + Pav (A-1) 


Po = Patm + Pep (A-2) 


el 


- Numbers in brackets refer to the numbered references listed at the end of 


this Documentation File. 
A4 
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During the thermal breathing process, the pressure, volume and temperature 
vary from a certain combination of values at the minimum condition to be certain 
combination of values at the maximum condition. At the minimum condition 1, it 
is assumed that all of the variables are simultaneously at their minimum values; 
and at the maximum condition 2, it is assumed that all of the variables are 
simultaneously at their maximum values. The value of the variables at the 
minimum condition 1 and maximum condition 2 are listed in Table Al. 


Table Al - Value of the Variables at the Minimum and Maximum Conditions 


Minimum Maximum 
Variable Condition 1 | Condition 2 
Gas space total pressure 
Atmospheric pressure 


Gas space gage pressure 


Stock vapor pressure 


Air partial pressure 


Gas volume 
Gas temperature 


Liquid surface temperature 


From the ideal gas law, the total number of moles of gas, ny, in an enclosed 
volume, V, at temperature, T, and pressure, P, is given by: 


=— (A-3) 


Assume that the gas mixture in the tank vapor space “is a two-component 
mixture consisting air and stock vapor. The mole fraction of air, ya, and the 
mole fraction of stock vapor, yy, may be determined by Eqs (A-4) through (A-7): 


AS 
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n 
Wyrm (A-4) 
sf 
P 
A 
y co (A-5) 
Ap 
n 
vy" a (A-6) 
nT 
P 
vy" - (A-7) 


The mole fraction of air, yA, may be expressed in terms of the mole fraction of 
vapor, yy, as follows: 


yaz l- YW (A-8) 


During the thermal breathing process of the tank vapor space, the number of 
moles of air, na, in the volumes, V] and V2, is assumed to remain the same. This 
assumption may be expressed as follows: 


MAl = NA2 (A-9) 
We may substitute Eq (A-4) into Eq (A-9) to yield: 

YA1 "T1 = YAa2 "T2 (A-10) 
We may substitute Eq (A-3) into Eq (A-10) to yield: 

Yal fr = Yao | (A-11) 


A6 
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Solving for Va and using Eq (A-8), we may write: 


ao [2] ae 
Yao) (Po) Uy 
ee p 
eG . 
: 1 - Yyo) bP Uy 


p 
eel 
Pry (Pa) {te 
Vv, =V a) | 
204 po} (P.) {7 
, Pye] Pd 
Po 
Pp, - P 
y. -y, | 2| | 2 (A-14) 
2 Tp, - yo} UT 
27 yd Wy 


Using Eq (A-14), the volume change due to thermal breathing, AV, may be 
determined as follows: 


AV = Vo - Vy (A-15) 


Pp, - Pp | , 
w=, aera [e aa (A-16) 
Po - Pyad Ty | 


We may substitute Eqs (A-1) and (A-2) into Eq (A-16) to yield: 


A7 
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Pry + Poy - P 
ews atm * Pay y Fe| of (A-17) 
7 Pam * Pap > Py) “a 
(Py - Py,) - (Pap - Pay) (T, - 1) 
Wey, f+ Stee 1 2) ate 
Pau + Pep - Pye 7 


It is a convenient to define the terms APy, APg and ATy as follows: 


APy = Py2 - Py (A-19) 
APg = Ppp - Pey (A-20) 
(A-21) 


ATy = T2 - T) = Ty2 - Tvl 


We may substitute Eqs (A-19), (A-20) and (A-21) into Eq (A-18) and expand the 
terms to yield: 


- APy = AP aT 
AV = V, f + y_ 6 | f + “al -1 (A-22) 


Pam + Pap.- Py2 


AT AP,, - AP (AP, - AP.) AT: 
Wey, 1+ A FA) - a} (a-23) 
Par * Ppp > Pyod Par * Pap > Pya 


Since the terms (ATy/Ty}) and (APy - APg)/(Patm + Pap - Pyo) are small, 
their product can be considered negligable. Thus, the product term in Eq (A-23) 
can-be neglected. £&q (A-23) then simplifies to the following: 


at AP, - AP 
av-v, (2 +—L—_42 — (A-24) 
Ty, Paty * Pap > Pye 
AS 


Copyright por American Petroleum Institute 
Wed Feb 09 12:30:48 2005 


API MPMS*¥19.1D 93 MM 0732290 0511453 217 


A3.0 VAPOR SPACE EXPANSION FACTOR 


The vapor space expansion factor, Kp, is defined as the ratio of the volume 
change, AV, to the initial volume, Vj, as follows: 


ees (A-25) 
Vy 


Substituting Eq (A-24) into (A-25) we obtain: 


(A-26) 


A4.0 SIMPLIFIED EQUATIONS FOR THE VAPOR SPACE EXPANSION FACTOR 


Eq (A-26) may be simplified for ease of calculation. Sections A4.1 through 
A4.4 present various simplifications that can be made. 


A4.1 Neglectinc the Term Ppp 


It should be noted that Ppp is small (about 0.03 psi) compared to PATM 
(about 14.7 psia) and can be neglected in the denomination of Eq (A-26) to 
yield: 


ATy  APy - AP 
gin Ge Be (A-27) 


Tya Pam > Pye 


Ke 


A4.2 Replacement of Ty) With Tra 


In the first term of Eq (A-27), the minimum vapor space temperature, Ty1> 
is close to the daily average liquid surface temperature, Ta, since both are 
absolute temperatures. Thus, for ease of calculation, we can replace Ty} with 
Tia in Eq (A-27) to yield: 


AQ 
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(A-28) 


A4.3 Replacement of Pyo With Pyas 


For low vapor pressure stocks, the stock vapor pressure, Py, is small 
compared to atmospheric pressure, Pam. Thus, we may replace the stock vapor 
pressure -at the minimum liquid surface temperature, Py2, with the stock vapor 
pressure at the daily average liquid surface temperature, Pya, in Eq (A-28) to 
yield: 


(A-29) 


Pam ~ Pya 


Eq (A-29) appears as Eq 4 in Ref. A7. 


A4.4 Use of a Simplified Equation for the Vapor Pressure Range, APy 


The vapor pressure of the stock may be determined from Eq (A-30), where the 
vapor pressure function constants A and B must be selected for the particular 
stock [see Tables 4 and 5 in Ref. A7]). 

Py = exp [A - (8/T1)] (A-30) 


We can determine the slope of the vapor space pressure function by taking 
its derivative with respect to the liquid surface temperature, TL, as follows: 


—_.—,* (A-31) 


Eq (A-31) can be written in differential form as follows: 


Ald 
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x AT, | (A-32) 


Eq (A-32) gives the vapor pressure range, APy, in terms of the liquid surface 
temperature range, AT;. 


For most hydrocarbon liquids, the liquid surface temperature range, AT, , 
may be related to the vapor temperature range, ATy, as follows (see Eq (F-16) 
in Section F): , 


— ATE = 0.50 ATy (A-33) 


Substituting Eq (A-33) into Eq (A-32), we obtain the following simplified 
equation for the vapor pressure range: 


0.50 B Py 
sa ne mee ae 


L 


(A-34) | 


Eq (A-34) may be substituted into Eq (A-29) to yield: 


AT 0.50 B){ P AP 
Ke = Ey 1+ Fees a (A-35) 
uf Tra J Paty > Py (Pata ~ Pya) 


A4.5 Neglecting the Term APg 


For most atmospheric storage tanks, the breather vent pressure and vacuum 
settings are typically +0.03 psi and -0.03 psi, respectively. Thus, the term 
APp/(Patm - Pya) is small (about 0.002 for low vapor pressure stock) compared 
to the term ATy/Tyq (about 0.040 for ATy = 20°F and Tiq = 520°R). For these 
cases, the last term in Eq (A-35) may be neglected to yield: 


All 
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A 0.50 8 P. 
Ths Tra J Pate > Pa 
In comparing Eq (A-36) to Eq (A-26), we can see that the calculation process 


is greatly simplified because it is not necessary to determine all of the 
variables, Thy, Tu2> TLA» Pv» Pv2 and Pya, but only the variables Tia and Pya.- 


AS5.0 CONCLUSION 


Equation (A-29) was selected for use in calculating the vapor space 
expansion factor, Ke, in API Publication 2518, Second Edition [A7]. This 
equation was derived from the ideal gas law and from the pressure, temperature, 
and volume conditions that exist in the vapor space of a fixed-roof tank 
containing a volatile liquid stock. Equation (A-29) was developed from a more 
complete expression, Eq (A-26), by incorporating several simplifications that 
make the calculations more user friendly, with little loss in accuracy. 


: This document is part of the API standards development process and is 
intended for use by API committees ONLY. It shall not be reproduced or 
circulated or quoted, in whole or in part, outside of the cognizant API 
committee(s) except with the written approval of API. This draft API standard 
will be formatted and edited prior to API publication. Copyright ©1990. 
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Pressure-Vacuum 
Vent 


Heat 


Figure Al - Schematic of the Tank Vapor Space Heating Process and the 
Resulting Volume Expansion Due to Thermal Breathing 


Al3 
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NOMENCLATURE 
SYMBOL DESCRIPTION UNITS 
A Constant in the vapor pressure equation psia 
AL Area of the stock liquid surface ft2 
a Defined by Eq (B-15) , dimensionless 
B Constant in the vapor pressure equation OR 
b _ Defined by Eq (B-16) dimensionless | 
Cy Average vapor concentration in the vented gas lbm/sft3 
Tank diameter ft 
E Evaporation loss lbm/day 
: Ec Evaporation loss calculated 1bm/day 
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! Hy Tank vapor space outage ft 
K Overall mass transfer coefficient between ( lbmole ) 
the liquid surface and the vented vapor 
ft2 hr mole frac. 
Ke Vapor space expansion factor dimensionless 
Ks Vented vapor saturation factor dimensionless 
My Stock vapor molecular weight lbm/Tbmole 
my Moles of stock vapor vented Tbmole 
PATM Atmospheric pressure psia 
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R Ideal gas law constant, (10.731) ft3 psia 
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Volume of gas vented during a single daily 
thermal breathing cycle 


Daily average stock vapor concentration in 
the vented vapor 


Daily average saturated stock vapor 
concentration 


Gas density 
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UNITS 


ft3 
mole fraction 


mole fraction 
1bmole/ft3 


This section of the Documentation File to API Publication 2518, Second 
Edition, contains the development of the vented vapor saturation factor, Ks. 


daily average stock vapor concentration in the vented gas, yy, to the stock vapor 
concentration, yy, in equilibrium with the stock liquid surface at the daily 
average liquid surface temperature. 


Section B2 presents the derivation of a theoretical equation for estimating 
the vented vapor saturation factor that is based on an analytical model of the 
daily thermal breathing process. 


APT MPMS*19.1) 93 MM 0732290 OS114be 2cT 
B1.0 INTRODUCTION 
The "Vented Vapor Saturation Factor", Ks, is defined as the ratio of the 
Section B3 presents the development of a correlation for estimating the 
vented vapor saturation factor that is based on test data. 
| 


B2.0 VENTED VAPOR SATURATION FACTOR MODEL 


B2.1 Model Description 


Figure 81 is a schematic of a fixed-roof tank that is partially filled with 
a volatile liquid stock and equipped with a pressure-vacuum vent. During the 
daily thermal breathing cycle, the gas mixture in the tank vapor space is 
initially heated from its minimum condition to its maximum condition (see Section 
A of this Documentation File for additional detail). Vapor is vented from the 
tank vapor space when the pressure increases to the pressure setting of the 
pressure-vacuum vent. As the gas mixture in the tank vapor space is cooled from 
its maximum condition back to its minimum condition, air is admitted to the tank 
vapor space when the pressure decreases to the vacuum setting of the pressure- 
vacuum vent. 


Evaporation of stock occurs from the liquid surface as the stock tries to 
saturate the air that was admitted to the tank vapor space. Test data indicates 
that there is a region at the top of the tank vapor space under the pressure- 
vacuum vent where there is a significant concentration gradient. 


BS 
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During the daily thermal breathing cycle, the moles of stock vapor that are 
expelled from the tank vapor space, my, may be estimated by Eq (8-1): 


my = pg AV yy (B-1) 


where the term AV represents the volume vented during a single thermal 
breathing cycle (see Section A2, Eq (A-24)). 


The vented gas contains less stock vapor per unit volume than the gas near 
the liquid surface. The daily average stock vapor concentration in the vented 
vapor is yy (unsaturated), and the daily average stock vapor concentration near 
the liquid surface is yy° (saturated). 


During the daily thermal breathing cycle, stock vapor evaporates and rises 
upward from the area near the liquid surface to replace the stock vapor lost as 
gas is vented from the tank vapor space. Stock will continue to evaporate as it 
tries to establish. a saturation condition at the top of the tank vapor space. 


The moles of stock that evaporate during a daily thermal breathing cycle may 
be estimated by Eq (B-2): 


my = KA, ty (Yy” - Yy) (B-2) 


where K is the overal] mass transfer coefficient between the liquid surface and 
the vented vapor. . 


After a series of repeated daily thermal breathing cycles where the same 
meteorological conditions occur, the stock vapor concentration in the vented 
vapor will vary during each thermal breathing cycle in a repeated manner, and the 
daily average stock vapor concentration in the vented vapor, yy, will achieve a 
steady value. This concentration value depends upon the rate at which the stock 
vapor lost from the tank vapor space is replaced by stock evaporated from the 
liquid surface. 
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The steady value of the daily average stock vapor concentration in the 
vented vapor, yy, may be determined by equating Eq (B-1) with £q {B-2) and 
solving for yy as follows: 


Pg AV = KA ty (yy - Yy) 


pp AV 
Vy 3 | = (yy? = Wy) (B-3) 


It is useful to define the saturation parameter, S, as follows: 


AV 
Sa ** | (B-4) 
KA, t 


Using this defined saturation parameter, Eq (B-3) may be written as: 


: vy $= (vy? - yy) (8-5) 
Eq (8-5) may now be solved for yy to yield: 
| | (B-6) 
Yy = ¥ 5 
Y V 1+5 


B2.2 Vented Vapor Saturation Factor Definition 


The vented vapor saturation factor, Ks, is defined by Eq (B-7) as the ratio 
of the daily average stock vapor concentration in the vented vapor, yy, to the 
daily average saturated stock vapor concentration, yy°. 


B7 
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k= ~ (B-7) 
fo] E 
dy 


When Ks = 1, the vented gas is completely saturated; when Ks = 0, the vented gas 
contains no stock vapor. 


82.3 Saturation Parameter 


ee 


The saturation parameter, S, as defined by Eq (B-4) may be written in terms 
of other evaporation loss parameters. First, note Eqs (B-8) threugh (B-11) as 


follows: . 
P 
oe am (8-8) 
R Ty 
2 
aD 
A, = | (B-9) 
AV = Vy Ke (see Eq (A-25) in Section A) (B-10) 
2 
DY H 
Vy = F y (B-11) 
. 4 | 


Substituting Eqs (B-8) through (B-11) into Eq (B-4), we obtain: 


p H 
$= [an | Ke (B-12) 
Ry) kt 
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The vapor space expansion factor, Kg, may be expressed in the following 
simplified form (see Eq (A-36) in Section A5): 


AT. 0.50 B P 
Ke = = 1+ lee (B-13) 
Tha Thy J Pam > Py | 


Substituting Eq (B-13) into Eq (B-12) and replacing Tyq with TLA (since both are 
absolute temperatures), we. obtain: : 


P H AT 0.50 B P 
a ccrlccs ad va Coe sai 
R Ta Kt Tia Tha Pam ~ Pya 


It is convenient to define the following dimensionless parameters a and b: 


Pary Hy AT 

— fore (B-15) 
Rt) K Tp 
0.50 8 

b= (B-16) 
Tha 


We also know that the daily average saturated stock vapor concentration, yy®, can 
be expressed as follows: 


P 
yo faa (B-17) 
p 
ATW 


B9 
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Substituting Eqs (B-15), (B-16) and (B-17) into Eq (B-14), we obtain: 


y fo] 
Sa |l+ | V 5 (B-18) 
] = Yy 


B2.4 Vented Vapor Saturation Factor Development 


Substituting yy from Eq (B-6) into Eq (B-7) we obtain: 


1 
K. = (B-19) 
s f + 1 


Equation (B-19) shows that as the saturation parameter, S, increases, the vented 
vapor saturation factor, Ks, decrease toward 0. Conversely, as S decreases, the 
value of Ks increases toward 1. 


Using Eq (B-18), the vented vapor saturation factor, Ks, may be written as: 


(B-20) 


Equation (B-20) shows that the vented vapor saturation factor, Ks, depends upon 
only 3 parameters: a, b and yy°. 


Inserting the expressions for a, b and yy® from Eqs (B-15), (8-16) and (B- 
17), we obtain the following final expression which contains ali of the 


variables: 
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(B-21) 


Po Hy AT 0.50 B 
ohare yr | 


2 
Rty K Tha Tha 


It should be noted that Ks will tend toward 1 as Hy tends toward 0. Also, Ks 
will tend toward 0 as Pyq tends toward Paty. 


Insufficient information is currently available to determine the overall 
mass transfer coefficient, K, and thus &q (B-21) was used only as a guide to show 
the dependancy of Ks on Pya, Hy and the other variables. 


Although it is possible to improve the above simplified analysis to develop 
a theoretical relation for the vented vapor saturation factor, Ks, it was decided 
instead to develop a correlation for Ks based on actual test data, as described 
in Section B3. However, the above simplified analysis was used as a guide in 
selecting the analytical form for the correlation equation and in selecting the 
parameters to include in the correlation. 


B3.0 VENTED VAPOR SATURATION FACTOR CORRELATION 


This section summarizes the development of a correlation for estimating the 
vented vapor saturation factor, Ks. 


Section B3.1 summarizes the saturation factors that were calculated from the 
API[38}*, EPA(20] and WOGA[17] test data. The API test data showed that the 
vented gas was near saturation conditions at all times. The EPA and WOGA Test 
Data, however, showed that the vented gas was not saturated, with the degree of 
saturation being less with increasing product vapor pressure, Pya, and increasing 


vapor space outage, Hy. 


= Numbers in brackets refer to the numbered references listed at the end of 


this Documentation File. 
Bll 
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The fact that the vented gas is not saturated with stock vapor reflects the 
effect of mass transfer rate limitations from the liquid surface to the area 
below the pressure-vacuum vent. As the vapor space outage increases, the 
distance that the stock vapor must travel from the liquid surface to the vent is 
lengthened. This lengthened distance decreases the mass transfer rate and thus 
the concentration in the vented gas. For high vapor pressure stocks, since the 
amount of stock vapor lost in each daily thermal breathing cycle is larger, a 
higher rate of evaporation from the liquid surface is required to replenish the 
stock vapor that is lost. Mass transfer rate limitations, however, limit the 
ability of the stock to replenish the vented vapor at these higher vapor 
pressures and thus reduce the degree of saturation in the vented gas. 


Section 83.2 presents the development of a correlation for the vented vapor 
saturation factor based on only the API and EPA test data. This correlation 
showed trends that are similar to those predicted by the theoretical analysis 
(see Eq (B-19)) in that the saturation factor approaches 1 as the vapor pressure 
or the outage approach 0, and the saturation factor becomes smal] as the vapor 
pressure or outage increase. 


Section 3.3 presents the development of a correlation for the vented vapor 
saturation factor based on the combined set of API, EPA and WOGA test data. This 
correlation showed the same trends as the correlation that was based on only the 
API and EPA test data, but there was more scatter of the WOGA test data from the 


correlation. 


B3.1 Summary of API, EPA and WOGA Test Data 


Table Bl summarizes the 10 API tests [38] along with the calculated 
saturation factor, Ks. The saturation factor for the AP] test data are very 
close to 1, with an average value for the 10 tests of 0.964. 


Table 82 summarizes the 15 EPA tests (20]. It was found that all of the EPA 
test data was suitable for use in calculating a saturation factor, with the 
exception of Tests EPA-4A, EPA-48 and EPA-4C. The reason for rejecting these 
tests is stated at the bottom of Table B2. 
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Table B3 summarizes the EPA test data along with the calculated saturated 
factor, Ks, for those tests which were selected in Table B2. Since the average 
liquid surface temperature was not measured during the EPA tests, the equation 
indicated in Note 5 at the bottom of Table B3 was used to estimate the average 
liquid surface temperature. 


Table B4 summarizes the 44 WOGA tests [17]. Out of the total of 44 tests, 
21 were found suitable to calculate a saturation factor. The reasons for 
rejecting the other tests is noted at the bottom of Table 84. 


Table BS summarizes the suitable WOGA test data and the calculated 
saturation factor, Ks. Only the crude oi] tests were used to calculate a 
saturation factor. The vapor pressure at the daily average liquid surface 
temperature was calculated utilizing the equations noted at the bottom of Table 
BS. No such relationships were available for the distillate and fuel oil 
products used in the WOGA test program. 


B3.2 Saturation Factor Correlation of API and EPA Test Data 


A 


In general, it was found that there was a higher quality in the API test 
data [38] and EPA test data [20] than in the WOGA test data [17]. The combined 
set of API and EPA test data were used to develop a saturation factor 
correlation. 


Figure 82 presents the correlation, where the saturation factor, Ks, was 
found to be related to the product of the vapor pressure, Pya, and the vapor 
space outage, Hy. The vapor pressure characteristics of the stock used in the 
EPA tests were readily known because they were single component stocks. 

83.3 Saturation Factor Correlation of API, EPA and WOGA Test Data 
The correlation of the API and EPA test data shown in Figure B2 is 


satisfactory, but is based only on test data from fuel oi] (API test data) and 
single component liquid stocks (EPA test data). 


B13 
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Although it was found that, in general, the API and EPA test data were of 
higher quality than the WOGA test data, it was decided to develop a single 
correlation which was fit to the combined set of API, EPA and WOGA test data. 
This combined data set includes 34 data points that extend up to a PyaHy value of 
about 78 psia ft. and include WOGA test data on crude oil. 


Table B6 summarizes the 34 data points which were used to develop the 
saturation factor correlation from the combined set of API, EPA and WOGA test 
data. 


Figure B3 is a plot of ((1/Ks) - 1) versus PyaHy for the API, EPA and WOGA 
test data. The test data were fit with a least squares linear correlation, as 
noted on Figure B3. The correlation coefficient, r2, was 0.76. Eq (B-22) is the 
resulting correlation. 

1 


aenweMnsee ees (B-22) 


Ks, = 
i + 0.0525 Pyally 


Figures B3 and B4 illustrate that the WOGA test data has more scatter in 
comparison to the API and EPA test data. Part of this scatter is believed to be 
due to the more uncertain vapor pressure characteristics of the stocks used in 
the WOGA tests and the fact that only a few vapor samples were taken during each 
WOGA test. 


Figure 84 illustrates the results of the correlation developed in Figure 83, 
where the saturation factor is plotted. versus Pyaly. The correlation indicates 
that the saturation factor approaches 1 as the vapor pressure or the vapor space 
outage approach 0. The correlation also shows that the saturation factor becomes 
small as the vapor pressure or the vapor space outage increase. These trends are 
in agreement with those that are predicted by Eq (B-21) of the theoretical 
analysis. 
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B4.0 CONCLUSION 


Equation (B-22) was selected for use in calculating the vented vapor 
saturation ‘factor, Ks, in API Publication 2518, Second Edition [A7]. This 
equation was developed from a correlation of the API, EPA and WOGA test data, and 
exhibits the same trends with varying Pya and Hyg that were exhibited by Eq (B- 
21) of the theoretical analysis. 
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Table B2 - Summary of EPA Tests [20] Selected 


Reasons 
for 
Rejection 


Selected 


Insulated 
Tank 
(Yes/No) 


Product Type 


Isopropano] 


Isopropanol 
Ethanol 
Ethanol 
Ethanol 


Glacial Acetic Acid 


Glacial Acetic Acid 


Formaldehyde 


Formaldehyde 


Formaldehyde 


Ethyl! Benzene 


Ethy] Benzene 


Cyclohexane 


= = © © ~ ~ « wes zeeZzene2az= 
<< «< «< @2 B@ Be mM wm we ae ut el 


Cyclohexane 


Cyclohexane 


Reasons for Rejection: (1) Tank is insulated. 
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Table B4 - Summary of WOGA Tests [17] Selected 


Product Insulated Vapor Selected Reasons 
Type Tank Outflow for 
(Yes/No) (sft3/day) | (Yes/No) | Rejection 


WOGA-1A 
WOGA-18 
WOGA-1C 
WOGA-2A 
WOGA-2B 
WOGA-3A 
WOGA-38 
WOGA-4A 
WOGA-4B 
WOGA-5 

WOGA-6 

WOGA-7A 
WOGA-7B 
WOGA-7C 
WOGA-8A 
WOGA-8B 
WOGA-8C 
WOGA-9A 
WOGA-9B 
WOGA-9C 
WOGA-10A 
WOGA- 108 
WOGA- 10C 
WOGA-11A 
WOGA-11B 
WOGA- 12 

WOGA-13A 
WOGA-13B 
WOGA- 14 

WOGA-15 


Crude Oi] 
Crude Oi] 
Crude Oi] 
Crude Oi] 
Crude Oi] 
Distillate 
Distillate 
Crude Oi] 
Crude Oi] 
Fuel Oi] 
Distillate 
Crude Oi] 
Crude Oi] 
Crude Oil 
Crude Oi] 
Crude Oi] 
Crude Oil 
Crude Oi] 
Crude Oi] 
Crude Oi] 
Distillate 
Distillate 
Distillate 
Crude Oi] 
Crude Oi] 
Crude Oi] 
Crude Oil 
Crude Oi] 
Distillate 
Crude Oi] 


. e e * 
WWWRP PD & 


OFT OFT OF bt ot be GC) 
o 


° 


ps 
Wi oaomoad0di 
. 2 @ 


ry 
@i Minne 


° 


. et et et ee ee ee ee et PP 2 2S 2 2S SS 


mw 
eOoOoonreooeo 


| il amd 


a - - 2. 2. 2 2. ? ? . 2) 2. 222 es 2 2S 2 eS 2 
Ne awe +s NNN UU 


Copyright por American Petroleum Institute 
Wed Feb 09 12:30:51 2005 


API MPMS*1L9-1D 93 MM 0732290 0511477 750 


Table B4 - Summary of WOGA Tests {17] Selected (Continued) 


Insulated Vapor Selected Reasons 
Tank Outflow for 
(Yes/No) (sft?/day) (Yes/No) | Rejection 


Crude Oi] §.2 Y 
Crude Oi] 5.2 Y 
Crude 0i1 §.2 Y 
Crude 0i1 5.5 Y 
Crude 0i1 5.5 Y 
Crude 0i1 0.4 N 
Crude 0i] 0.4 N 
Crude 0i1 0.4 N 
Crude Oi] 3.0 Y 
Crude 0i1 3.0 Y 
Crude 011 4.1 Y 
Crude 0i1 4.1 Y 
Distillate --- 


Distillate 


‘ 
‘ 
’ 


WOGA-21B 


Reasons for Rejection: 


(1) Tank is insulated. 
(2) Product is a distillate or fuel oi]. The vapor pressure versus 
temperature behavior of these stock types was not well established 


in these tests. 
(3) The liquid sample may not be representative of the tank stock during 


the test. 
(4) The RVP was less than 1.0 psi, and the stock vapor pressure is, there- 


therefore, questionable. 
(5) The vented gas volume outflow, Q, appears to be in error. 
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Table B6 - Summary of Test Data Used to Develop the 
Vented Vapor Saturation Factor Correlation 
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Vented Vapor 


eanwnd 


i Bl - Schematic of the Tank Vapor Space for the 
ail Vented Vapor Saturation Factor Analysis 
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i B2 - Saturation Factor, Ks, Versus PyaHy 
oi for the API and EPA Test Data 
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Figure BS - Saturation Factor, Ks, Correlation 


A EPA Test Data 
QIWOGA Test Oata 


© API Test Data 
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NOMENCLATURE 


DESCRIPTION 


Area 
Defined by Eq (C-35) 


Ambient temperature range factor, defined by 
Eq (C-44) 


Heat capacity 


_ Diameter 


Vapor space aspect ratio, defined by Eq (C-77) 
Defined by Eqs (C-19), (C-20) and (C-21) 

Solar insolation factor, defined by Eq (C-78) 
Defined by Eqs (C-22), (C-23) and (C-24) 


Heat transfer coefficient ratio, defined 
by Eq (C-79) 


Defined by Eqs (C-25), (C-26) and (C-27) 


Height 


Daily total solar insolation on a horizontal 
surface 


Heat transfer coefficient 

Solar insolation intensity 

Solar absorptance ratio, defined by Eq (C-80) 
Defined by Eq (C-68) 


Vapor space temperature factor, defined by 
Eq (C-60) 


Mass 

Heat transfer rate 

Solar reflectivity of the ground 

Ratio of Iq to Hy, defined by Eq (C-71) 
Temperature , 

Temperature change 

Time 

Overall heat transfer coefficient 


Zenith angle, or the angle between the sun and 
the normal to a horizontal surface 


Solar absorptance 
Atmospheric solar transmittance 
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UNITS 


ft2 
B/hr ft2 


dimensionless 
B/1bm OF 

ft ; 
dimensionless 
dimensionless 
dimensionless 
dimensionless 


dimensionless 
OF 
ft 


B/day ft2 
B/hr ft? OF 
B/hr ft2 
dimensionless 
dimensionless 


dimensionless 
Tbm 

B/hr 
dimensionless 
day/hr 

OF 

OF 

hr 

B/hr ft2 OF 


deg. 
dimensionless 
dimensionless 
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SUBSCRIPTS 
A Ambient (or air) 
AV Average 
B Beam component of solar insolation 
D Diffuse component of solar insolation 
G Gas 
GD Ground 
H Horizontal surface solar insolation 
I Inside 
L Liquid 
MN Minimum 
MX Maximum 
N North shell (or the half of the tank shell 
facing away from the sun) 
NI North shell inside 
NO North shel] outside: 
0 Outside 
_R Roof 
RI Roof inside 
RO Roof outside 
s) South shell (or the half of the tank shell 
facing toward the sun) 
SI Soutn shell inside 
SN Solar noon 
$0 South shell outside 
T Temperature 
| 
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C1.0 INTRODUCTION 


This section of the Documentation File to API Publication 2518, Second 
Edition, presents the development of the equation for the vapor space temperature 
factor, Ky, which is defined as the ratio of the gas space temperature range, 
ATg, to the ambient temperature range, ATa- 


Section C2 describes the heat transfer model. The major assumptions used in 
the mode] development are: 


0 


Assumption 1 
The gas space is fully mixed (i.e. it is at a uniform temperature and 


composition). 


Assumption 2 . 
The liquid space is fully mixed (i.e. it is at a uniform temperature 


and composition), and remains at a constant temperature during the 
daily cycle. 


Assumption 3 
The tank wall in the gas space can be treated as three (3) separate 


elements: (1) the roof; (2) the half of the tank wall facing away from 
the sun; and (3) the half of the tank wall facing the sun. Each tank 
wall element can be characterized by a single temperature, which varies 
during the daily cycle. 


Assumption 4 
The affects of rain and snow precipitation are not included in the 


model. 


These assumptions are the same as those upon which the API Computer Model is 
based [30, 38]”. 


LS 


- Numbers in brackets refer to the numbered references listed at the end of 


this Documentation File. 


C6 
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Section C3 presents the heat balance differential equations for each of the 

tank wall elements and the gas space (see Eqs (C-1),. (C-3), (C-5) and (C-7)). 

- These ordinary differential equations are essentially the same as those used in 

the API Computer Mode? [30, 38], where they are solved by stepwise numerical 
integration. 


A sensitivity analysis [30] of the API Computer Model showed that the gas 
and tank wall heat capacity terms in the differential equations have a negligible 
affect on the computed results. Thus, the following additional assumption was 
made: 


1) Assumption 5 
The heat capacity terms in the energy _balance equations can be 


neglected in comparison to the other heat transfer terms. 


With this assumption, the differential equations reduce to a set of four (4) 
simultaneous algebraic equations (see Eqs (C-16), (C-17), (C-18) and (C-29)). 


Section C4 presents the solution to these simultaneous equations by. solving 
for the gas temperature, Tg (see Eq (C-36)), and gas temperature range, ATg (see 
Eqs (C-58) and (C-59}). 


Section C5 presents the vapor space temperature factor equation for the 
general case of a two color fixed-roof tank, where the roof and shell -are painted 
different colors (see Eqs (C-74) through (C-80)). 


Section C6 presents the vapor space temperature factor equation for the 
simplified case of a two color fixed-roof tank where typical solar insolation 
parameters are used (see Eqs (C-89) through (C-94)). 


Sections C7, C8 and C9 present the vapor space temperature factor equation 
for progressively simplified cases where more typical information is used instead 
of detail information for a particular tank (see Eqs (C-97), (C-106), (C-110) and 
(C-115)). As less detail information is used for the calculation of Ky, the 
estimation equations become simpler, but the estimated value of Ky becomes less 


accurate for a specific tank. 


C7 
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C2.0 HEAT TRANSFER MODEL DESCRIPTION 


Figure Cl is a schematic of the energy flows and temperatures for a fixed- 
roof tank. The gas space tank wall is divided into three (3) tank wall 


elements: 


) the roof, 

0 the half of the gas space tank shell that faces away from the sun 
(referred to herein as the “north shell"), and 

o the half of the gas space tank shell that faces toward the sun 
(referred to herein as the “south shell"). 


Each of these tank wall clements is characterized by a single, different 
temperature, Tp, TN and Ts, respectively. The temperature of each wall element 
varies with time over the course of the daily cycle in response to a heat balance 
on the wall element. 


The elements exchange heat on both their inside and outside surfaces. The 
inside of each element exchanges heat with the gas, which is characterized by a 
single temperature, Tg, by natural convention heat transfer. Heat transfer on 
the inside surface by long wave length thermal radiation, however, is neglected 
because the magnitude of this heat transfer rate is small in comparison to the 
natural convection heat transfer rate. 


The outside of each element exchanges heat with the ambient air by 
convection due to the wind. In addition, the outside of each element receives 
solar insolation due to beam solar insolation (except the north shell), diffuse 
solar insolation, and ground reflected solar: insolation. Each wall element also 
exchanges heat by long wave length thermal radiation with the tank surroundings. 
This last mode of heat transfer is incorporated in the outside heat transfer 
coefficient for each element, hpyg, hyo and hso, respectively. 


It should be noted that the thermal resistance of any thermal insulation 
material applied to the outside surface of the roof and shel] elements can be 
incorporated in the outside heat transfer coefficients hpo, hyo and hso.- 


cs 
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The gas space exchanges heat by natural convection with the wall elements 
and the liquid surface, which is assumed to be at a constant liquid temperature, 
Ti, during the daily cycle. The gas space is assumed to be fully mixed, so that 
it can be characterized by a single temperature, Tg. 


C3.0 HEAT BALANCES 


Sections C3.1 through C3.4 describe the heat balances that can be written 
for each of the four elements. 


C3.1 Roof Heat Balance 


(C-1) 


The term on the right hand side of Eq (C-1), which represents the heat capacity 
effect of the roof, will be neglected because its magnitude at most times during 
the course of the daily heating cycle is small in comparison to the other terms. 
Eq.(C-1) then reduces to: 


gro + Gri = 9 (C-2) 
C3.2 North Shell Heat Balance 
. dTy 
Gno * Sut * "NoN ne (¢-3) 


The term on the right hand side of Eq (C-3) will again be neglected because its 
magnitude is small in comparison to the other terms. €q (C-3) then reduces to: 


ano + GNI = 0 (C-4) 


cg 
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C3.3 South Shell Heat Balance 


Oe 


dT, ie 
Aso * Igy = MSs ra (C-5) 


The term on the right hand side of Eq (C-5) will again be neglected because its 
magnitude is small in comparison to the other terms. fq (C-5) then reduces to: 


aso + ast = 0 (C-6) 
C3.4 Gas Space Heat Balance 
dT 
Qpy + Wr t 857 +H ™ “MEG [OL (C-7) 


The term on the right hand side of Eq (C-7) will again be neglected because its 
magnitude is smal] in comparison to the other terms. Eq (C-7) then reduces to: 


ari + GNI + ast + GL = 9 (C-8) 


¢3.5 Heat Transfer Equations 


Eqs (C-9) through (C-15) are the heat transfer equations for each of the 
terms RQ) GRI> QNO» GNI» SO» 4SI and q_ that are used in Eqs (C-1) through 


(C-8). 
“aRO = heoAe (Ta - Tr) + an (Anip Cos @ + ArIp) (C-9) 
ari = brrAr (Te - TR) {C-10) 
1 1 
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ant = SytAn (Tg - TN) (C-12) 


] 
= hecAe (T. - Te) + Qe |DHnIp Sin @ + — Acly, + 
Asg * NggAs (Ty - Tg) + ag JOMGE > stb 
3 Acl'an (Ip Cos @ + I,) (C-13) 
> SYeD ‘'B ny 
ast = hsjAs (Te - Ts) (C-14) 
au = hLAL (Te - 11) (C-15) 


C3.6 Tank Element Temperature Equations 


Eqs (C-9) through (C-15) can be substituted into the tank element energy 
balance relations, Eqs (C-2), (C-4), and (C-6), to result in the following tank 
element temperature equations, Eqs (C-16) through (C-18). 


Tr = epTg + fala + OR . (C-16) 
Ty = enlTg + fyTa + ON (C-17) 
Ts = esTg + fsTa + 9S (C-18) 


where ep, en, €s, fr, fn, fs, GR» ON and Qs are defined by Eqs (C-19) through 
(C-27). 


@ = —— (C-19) 
hey + hp, 


Cll 
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| het : (C-20) 
ey * Se ; 
Dwr + hy 
; -| hs | (c-21) 
S 

Ney + hg 
f -| Ro ] (c-22) 
R 

hey + hp 
Pe bno | (C-23) 
N 

Dar + by 
f. = [= (C-24) 
S 

hey + he 


@p (Ip Cos @ + Ip) (c-25) 
ae bea 


hay + MR 


i, 
ay [> ty + 5 Teo (lg os @ * Ip) 


-26 
3, -|—————______—— (C-26) 


hwr * xo 
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2 1 1 
a. |- I, Sin @+—- 1) + - Ven (I Cos @ + I,) 
s) ; B 2 D2 GD ‘'B D 
95 * ee (C-27) 


(4.0 VAPOR SPACE TEMPERATURE FACTOR 


C4.1 Gas Space Temperature Equation 


The gas space temperature equation can be developed from the gas space 
energy balance equation, fq (C-8). Substituting Eqs (C-10), (C-12), (C-14) and 
(C-15) into Eq (C-8), we obtain: 


hrrAr (Tg - Tr) + bwtAn (Te - TN) + 
hsjAs (Tg - Ts) + HLAL (Te - TL) = 0 (C-28) 


Eq (C-28) may be rewritten as follows: 


Tg (hprAr + hyrAN + hszAs + HALAL) - 
Tp (hprAr) - Tx (hwrAN) - 
Ts (hsqAs) - TL (hLAL) = 9 (C-29) 


Eqs (C-16), (C-17) and (C-18) may now be substituted into Eq (C-29) for the terms 
Tr. Ty and Ts, respectively. The resulting equation may be rearranged by 
factoring out the temperatures Tg, Ta and Ty to obtain: 


Tg [hgrAp (1 - ep) + BNIAN (1 - eN) + hsyAg (1 - es) + HLAL) = 
Ta (hRTARFR + hNIANfn + HsTAsfs) + TL (hyAc) + 
(hetArgr + hyTANSN + hstAsgs) (C-30) 


Substituting the expressions for eR, eN, eS: fp, fy. fs, OR» SN and 9S found in 
Eqs (C-19) through (C-27), the above Eq (C-30) becomes: 
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Tg (UpAR + UNAN + UsAs + ALAL) = 


Ta (UpAR + UNAN + UsAs) + Tr (ALAL) + BAL (C-31) 
where Up, Uy, Us and B are defined by Eqs (C-32) through (C-35). 
ho,h 
4 [ee (C-32) 
hey + hp 
hy ph 
> fe | (c-33) 
Ayr + hy 
he, h 
vs [rae (C-34) 
hgy + hg 
A 
B = alr Gp (Ip Cos 6 + Ip ) + 
h A 
R L 
n| [An gi aoe Geet ceoanrats) 
—j} I—| ay I~ Ip + - I, Cos @ + + 
h A, N Io dD’ 3 GD ‘B 9) 
Us Ay 2 1 1 
| |] es 7 Ip Sin@+- Ip + ‘ Tep (Ip Cos @ + Ip) (C-35) 
h A a 2 2 
SY SL 
Eq (C-31) may now be solved for Tg to obtain: 
T,(UpA, + UyAy + UcAc) + T, (h, A, ) + BA 
| Tatra * nin * sis) bk L (c-36) 


G (UpAp + UnAy + UcAe + hy Ay ) 
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C4.2 Maximum Gas Space Temperature 


The maximum gas space temperature, 1G, MX» will generally occur at Solar 
Noon. &q (C-36) then becomes: 


Ta sw(UpAn * Un + UsAs) + TUMLAL) + BswA | (C-37) 


T 
G, WX (Upha + Uyay + Ushs + bpAL) 


where Boy represents Eq (C-35) with the terms Ig, Ip and @ evaluated at Solar 
Noon. 


C4.3 Minimum Gas Space Temperature 


The minimum gas space temperature, Tg my. will occur in the night at the 
time of minimum ambient temperature, TAa,MN> when there is no solar insolation 
(i.e. when B= 0). Eq (C-36) then becomes: 


_ Tanne + Unay + UsAs) + TH (MAL) (6-38) 


T 
G,MN (UpAe + Uvby + Ushs + WAL) 


C4:4 Average Gas Space Temperature 


The average gas space temperature, Tg,ay» is defined by Eq (C-39): 


T +T 
G,MX G,MN 
2. —L. Le = 
Te VAY ; (C-39) 


Substituting Eqs (C-37) and (C-38) into Eq (C-39), we obtain: 


1 
Ta ay (peg + UyAy + UsAs) + TLChLAL) + psn 


(UpAp + UyAN + UsAs + hy Ay) 


G,AV (C-40) 


C15 
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C4.5 Gas Space Temperature Change 


The gas space temperature change, ATg, may be determined by subtracting 
Eq (C-38) from Eq (C-37) to yield: 


ATg = (Te,mx - TG,MN) (c-41) 
ATy sw(Upap + UnAw + UsAs) + Boni : 
(UpAp + UnAy + Ugg + AL) 
where ATa sn is are defined as follows: 
(C-43) 


ATa,sn ® (Ta,sN - TA,MN) 


C4.6 Ambient Temperature Range Factor 


It is convenient to define the ambient temperature range factor, C, as 


follows: 
AT 
Ce —ALSN (C-44) 
AT, 
or 
ATa,sn = C ATa (C-45) 


where ATa is defined as follows: 


ATa = (Ta,mx - TA,MN) (C-46) 


Figure C2 depicts a sinusoidally varying daily ambient temperature. The 
ambient temperature varies from its minimum value, Ta,MN» at time tn to its 
maximum value, Ta Mx, at time tyx. This sinusoidal variation is described by Eq 


(C-47). 


C16 
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AT, 360 
Ty = Ta AV + i Sin r¢ (t - tuy + 6) (C-47) 


It is common for the maximum ambient temperature to occur about 2 hr after Solar 
Noon, or at 14:00 hrs. For this reason, we will select: 
tuy = 14 hr (or 14:00 hrs) (C-48) 


With this value for twx, Eq (C-47) becomes: 


AT, 360 ( , (c-49) 
Tel + — Sin ||— | (t - 8 - 
Solar Noon occurs at tsy: 


tsy = 12 hr (or 12:00 hrs) (C-50) 


The ambient temperature at Solar Noon, Ta,sn, from Eq (C-49) is: 


: 7 AT, 360 ( 
= + —— Sin ||—j (12 - 8) 
A,SN A, AV 2 24 


Ty ay + 0-433013 AT, (C-51) 


T 


A,SN ” 'A,AV 


we can now calculate C from Eq (C-44) using Eqs (C-43) and (C-51) as follows: 


AT 
c= -ALSN 
| AT, 


C17 
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c = JASN ~_"A,IN 
AT, 


eee Ty av + 0.4330134T, - Th MN 


AT, 


AT, 


C = 0.933013 (C-52) 
C4.7 Area Equations 


For a fixed-roof tank with a flat roof, the areas Ap, An, As and AL may be 
determined from Eqs (C-53) through (C-56). 


id | (C-53) 
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C4.8 Gas Space Temperature Change Equation 


Using the area relationships of Eqs (C-53) through (C-56) and Eq (C-45), we 
may rewrite Eqs (C-40), (C-42) and (C-35), respectively as follows: 


|e (C-57) 


2H. 2He ‘ 
U, + Uyj—] + Uel[-—] + 
RON 5 S] 5 L 


(C-58) 


Bon = | Qp (Ip Cos 8 + Ip) + 


Uy 2H. l 1 a 1) 
j——}] [—] a |- I, + - 7 Cos 6 + + 
hy D N 2 D 2 GD ‘'B D 

Pa ae ee 


SN 
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The SN at the end of Eq (C-59) indicates that the terms Ig, Ip and @ are 


evaluated at Solar Noon. 


C4.9 Vapor Space Temperature Factor 


The vapor space temperature factor, Kt, is defined. as the ratio of ATg to 
ATa, as follows: 


AT, 
Ky = |— (C-60) 
AT 


Using Eq (C-58), K7 becomes: 


2H 2H B 
D D 


_—————— : (C-61) 


2H 2H 
D 1) 


where Bsy is determined from &q (C-59). 
C5.0 TWO COLOR TANK 


If the tank roof is painted a different color than the tank shell, then the 
solar absorptance of the roof is ap and the solar absorptance of the shell is 
as, where: 


aN = as (C-62) 
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To simplify Eqs (C-57), (C-59) and (C-61), we will assume that: 


C-63 

ho = hro = bno = ‘so (C-63) 
C-64 

hy = heyy = buy = hsy = hy ( ) 
-65 

U = Up = UW = Us (C-65) 


Using these simplifications, Eqs (C-57) and (C-61) reduce to: 


(C-66) 
Ss 
He) oh 
2 +f} «3 
0 0 
(C-67) 
Secs coke See ip rns 
Tg AV 
h 
2 + 5 pt 
D hy 
where: 
ij Bsn (C-68) 
a UAT, 


C21 
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@ a H 
Key = for Fs ‘ rea (Ip Cos @ + Ip) + 
hoAT a> D 


| E I, Sing +! | (C-69) 
oilx 8 D 
SN 


It is convenient to rewrite Eq (C-69) in terms of Iy and Ip instead of Ig 
and Ip. Note the following relation, Eq (C-70), between Iy, Ig and Ip, and the 
relation, Eq (C-72), between Iq and Hy. 


Iy = Ip Cos 8 + Ip (C-70) 
I 
ae (C-71) 
Hoon 
H 
Iq = rH Hy (C-72) 


Using Eqs (C-70) and (C-72), Eq (C-69) may be rewritten as follows: 


2 I I . 
[5 | f : 4 Tan 9 + 2 (C-73) 
o| {lx ly ly | 


SN 


It ig convenient to define the dimensionless parameters E, F, G and J as 
described by Eqs (C-77), (C-78), (C-79) and (C-80), respectively. Using these 
defined parameters, Eqs (C-66), (C-67) and (C-73) may be rewritten as follows for 
the case of a two color fixed-roof tank. 
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Two Color Fixed-Roof Tank 


C(1 + 4E) + K 
eee eee eee (C-74) 


K 
T (2 + 4E + 6) 


1 
Tay (1 + 46) + T (1 + @) + ksh AT, 


a (C-75) 


T = 
G,AV (2 + 4€ + G) 


2 Ip Ip 
Key = F 4(d + 2rppE) + 2E{|-| [1 - —] Tan @ + — (C-76) 
SN GD . I I 
HI ISN 


where: 


(C-77) 
(C-78) 
(C-79) 


(C-80) 


C6.0 TYPICAL SOLAR INSOLATION PARAMETERS 


The expression for Ksy, Eq (C-76), may be simplified for the case of the 
typical solar insolation parameters used to generate the API Computer Data Base 
[39]. These typical solar insolation parameters are developed in Section D of 
this Documentation File and are based on the following assumed conditions: 
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Ground reflectivity, rep, of 0.1, 
Tank site latitude is 35 deg. north, 

_ Tank site elevation is 2000 ft. above sea level, and 
‘The year day number is 105 (April 15). 


o o ¢6 6 


For the entire API Computer Data Base, parameters 7), 7B. YH and @ are the 
following constant parameters (see Table D2): 


Tp = 0.0767938 (C-81) 
7p = 0.660790 (C-82) 
ry = 0.133277 (C-83) 
@ = 25.5848 deg. (C-84) 


The ratio Ip/Iy can be expressed in terms of 7p and 7g as follows: 


I T 
JO. Wy (C-85) 
ly Ta+T . 


Using the values of rp and 78, listed by Eqs (C-81) and (C-82), we can 
calculate the ratio Ip/Iy from Eq (C-85) as follows: 


Ip 0.0767938 
0.660790 + 0.0767938 


= 0.104115 (C-86) 


Substituting the above values for Ip/Iy, rep, TH and @ into Eqs (C-76) and 
(C-Ta), we obtain: 
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Ken = F (J + 0.954378 E) (C-87) 
0.133277 a. H 
F = ay (C-88) 


Eqs (C-52) and (C-87) may be substituted into Eqs (C-74) and (C-75) to yield the 
following result of Ky and Tg ay for the case of a two color fixed-roof tank with 
typical solar insolation parameters: 


Two Color Fixed-Roof Tank With 


Typical Solar Insolation Parameters: 


0.933(1 + 4£E) + F (J + 0.954E) 


Sa e— ee (C-89) 
(2 + 4E + &) 
1 
Ta AV (1 + 4E) +7, (1 +6) + shia (J + 0.954E) 
Te av * (C-90) 
ane (2 + 4& + G) 

where: 
(C-91) 
(C-92) 
(C-93) 
(C-94) 
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C7.0 ONE COLOR TANK 


If the roof and shell of the fixed-roof tank are painted the same color, the 
solar absorptance of the roof and shell paint can be represented by the same 


value, a: 

@ = as = aR (C-95) 
Substituting these values into Eq (C-94), we see that: 

J=] (C-96) 


Eqs (C-95) and (C-96) may be substituted into Eqs (C-89) and (C-90) to give the 
following result of Kt and Tg ay for the case of a one color fixed-roof tank with 
typcial solar insolation parameters: 


One Color Fixed-Roof Tank With 
Typical Solar Insolation Parameters: 


0.933(1 + 4€) + F(1 + 0.954E) 
(2 + 4E + G) 


1 
Ty ay (i + 4) + T, (1 + 6) + sat (1 + 0.954E) 


T = 
G,AV (2 + 4E + G) 


where: 
H 

ee ee 
8) 


Se - a “j 
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Thus, the vapor space temperature factor, KT, depends only on the three 
dimensionless parameters, £, F and G. 


Figures C3 through C6 illustrate the dependency of Ky on E, F and G. 
C8.0 TYPICAL HEAT TRANSFER COEFFICIENTS AND TYPICAL VAPOR SPACE ASPECT RATIO 
Figures C7 and C8 show the effect of heat transfer coefficient ratio, G, on 
the average percent difference of estimated emissions compared to the API 


Computer Data Base [39] and the API Test Data [33]. These figures show that the 
following average value may be selected for G: , 


G = 0.45 (C-102) 


This value of G corresponds to the following average values for the inside and 
outside heat transfer coefficients, hy and ho, respectively: 


hy = 0.65 B/hr ft OF (G-103) 

ho = 1.45 B/hr ft2 OF (C-104) 

Figure C5, which is for a value of G = 0.50, shows that Ky depends little on 
the value of the vapor space aspect ratio, E, as the insolation factor, F, varies 
from 1 to 10. An average value of E may be selected as follows: 

—E = 1.0 | (C-105) 

Substituting the typical values of G, hy, ho, and E listed in Eqs (C-102) 


through (C-105) into Eqs (C-89) through (C-94) for a two color fixed-roof tank, 
we obtain the following results for Ky, ATy and Tg,av: 
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Two Color Fixed-Roof Tank With 
Typical Solar Insolation Parameters, 
Typical Heat Transfer Coefficients, and 


Typical Vapor Space Aspect Ratio: 


aly, 
Ky = 0.723 + 0.0279|-—— (C-106) 
AT 
ATy = 0.723A4Ta + 0.0279aHy (C-107) 
Te, ay = 0.775TA,AV + 0.225T_ + 0.0140aHy (C-108) 
where: 
(C-109) 


ho and E listed in Eqs (C-102) 
lor fixed-roof tank, 


Substituting the typical values of G, hy, 
through (C-105) into Eqs (C-97) through (C-101) for a one co 
we obtain the following results for KT, ATy and Tg, av: 


One Color Fixed-Roof Tank With 
Typical Solar Insolation Parameters, 
Typical Heat Transfer Coefficients, and 


Typical Vapor Space Aspect Ratio: 


aH, 
= 0.723 + 0.0279 |= 


Al, 


Ky 


ATy = 0.723 ATq + 0.0279 chy 


Tg, AV = 0.775Ta ay + 0.225T, + 0.0140aH} 
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C9.0 TYPICAL METEOROLOGICAL PARAMETERS AND TYPICAL PAINT SOLAR ABSORPTANCE 


Figure C9 is a plot of the ambient temperature range, Ala, as a function 
of the daily. total solar insolation on a horizontal surface, Hy, for the selected 
U.S. locations listed in Table Cl [the meteorological data are from Table 6 in 
Ref. A7]. Figure C10 is a plot of the ratio (Hy/ATa) as a function of ATp for 
the same locations presented in Figure C9. For these locations, the average 
annual value of (Hy/ATa) is: 


H 
HH . 66.0 B/ft* day °F (C-113) 
AT, 


If specific information is not available on the tank paint color and paint 
condition, a white shell and white roof, with the paint in good condition, may be 
assumed. For this assumed condition, the following value of paint solar 
absorptance may be used: 


@ = 0.17 (C-114) 


Substituting the typical values for (Hy/4Ta) and @ from Eqs (C-113) and 
(C-114) into Eqs (C-110) through (C-112), we obtain the following typical result 
for Ky, ATy and Tg av. 


One Color Fixed-Roof Tank With 

Typical Solar Insolation Parameters, 
Typical Heat Transfer Coefficients, 
Typical Vapor Space Aspect Ratio, 
Typical Meteorological Conditions, and 


Typical Paint Solar Absorptance: 


Kt = 1.04 (C-115) 
ATy = 1.04 Ala (C-116) 


Te,av = 0.775Ta Ay + 0.2257, + 0.157ATA (C-117) 
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C10.0 CONCLUSION 


In summary, Equations (C-61), (C-74), (C-89), (C-97), (C-106), (C-110) or 
(C-115) may be used to estimate Kr, depending upon the level of detail in the 
information available for a specific tank. As less information is available, the 
equation used for the calculation of KT becomes simpler, but the calculated value 
of Ky will be less accurate for the specific tank. 


Equation (C-107) was selected for use in calculating the daily vapor 
temperature range, ATy, in API Publication 2518, Second Edition [A7]}. This 
equation is based on a comprehensive analytical heat -transfer model of the tank 
vapor space during a daily heating cycle. Equation (C-107) was developed from a 
more complete expression, Eq (C-58), by incorporating several simplifications 
that make the calculations more user friendly, with little loss in accuracy. 
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Table Cl - Annual Average Meteorological Data for Selected 
U.S. Locations [22,35] 


Location 


Birmingham, AL 
Montgomery, AL 


Santa Maria, CA 
Denver, CO 

Grand Junction, CO 
Wilmington, DE 
Atlanta, GA 
Savannah, 6A 
Honolulu, HI 
Chicago, IL 
Springfield, IL 
Indianapolis, IN 
Wichita, KS 
Louisville, KY . 
Baton Rouge, LA 
Lake Charleston, LA 


Grand Rapids, MI 
Minneapolis, MN 
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Table Cl - Annual Average Meteorological Data for Selected 
U.S. Locations [22,35] (Continued) 


(x) 
Ta,MX | TA,MN Hy Ta, av | ATA ATA 
Location (or) | (OF) |(B/day ft2)| (°F) | (OF) |(B/day ft? °F) 


Jackson, MS 
Billingsw, MT 
Las Vegas, NV 


Philadelphia, PA 
Pittsburg, PA 
Providence, RI 
Columbia, SC 
Sioux Falls, SD 
Memphis, TN 


Hungtinton, WV 
Cheyenne, WY 
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Figure Cl - Schematic of Energy Flows and Temperatures 
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Figure C2 - Sinusoidal Daily Ambient Temperature Variation 
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Solar Insolation, Hy (Btu/day ft2) 


Figure C9 - ATq Versus Hy for Selected U.S. Locations 
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API PUBLICATION 2518 
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SECTION D 
NOMENCLATURE 
SYMBOL DESCRIPTION UNITS 
A Altitude ft 
a Defined by Eq (D-12) --- 
ao Defined by Eq (D-7) oo 
ag” Defined by Eq (D-4) --- 
a) Defined by Eq (D-8) oe 
ay Defined by Eq (D-5) --° 
b Defined by Eq (D-13) --- 
Hy Daily total solar insolation on a horizontal surface B/day ft2 
IB, SN- Beam solar insolation at solar noon Bar ft? 
Ip, SN Diffuse solar insolation at solar noon B/hr ft2 
TH. SN Hourly total solar insolation on a horizontal 
, surface at solar noon B/hr ft2 

k Defined by Eq (D-9) mee 
k* Defined by Eq (D-6) ose 
n Year day number --- 
YO Correction factor for climate type from Table Dl --- 
Y] Correction factor for climate type from Table DJ --- 
YK Correction factor for climate type from Table 01 --- 
TH,SN - Ratio of Iy,sn to Hy --- 
é Declination angle deg. 
8SN Zenith angle at solar noon deg. 
7B,SN Atmospheric transmittance for beam solar insolation 

at solar noon --- 
7D,SN Atmospheric transmittance for diffuse solar insolation 

at solar noon --- 
¢ Latitude deg. 
w Sunset hour angle deg. 

Subscripts: 
B ' Beam 
8) Diffuse 
H Horizontal} 
SN Solar Noon 
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D1.0 INTROOUCTION 


This section of the Documentation File to API Publication 2518, Second 
Edition, contains equations that may be used to determine the solar insolation 
parameters required to calculate the vapor space temperature factor, Ky. These 
equations have been selected from Duffie and Beckman [24]* and are arranged in 
the order required to calculate Iy,sn» Ip,sn» and Ig,sn given n, ¢, Hy and A. 


Section D2 presents the equations to calculate each of the required solar 
insolation parameters. 


Section 03 presents sample calculations illustrating the use of the 
equations for a sample problem. Table D2 summarizes the results of the sample 
problem. These results are also used in Section C6 of this Documentation File to 
develop an equation for the vapor space temperature factor, Ky. 


D2.0 SOLAR INSOLATION EQUATIONS 


02.1 Declination, 6 


The declination, 5, is the angular position of the sun at solar noon with 
respect to the plane of the equator, where a north declinatiun is considered 
positive. The declination varies between -23.459 and +23.459 and is only a 
function of the day of the year, n, where 1 <n < 365. 


360 (284 + n) 
§ = 23.45 Sin |-————————_ (D-1) 
365 


[see Ref. 24, pg. 11, Eq (1.6.1)] 


eel 


* Numbers in brackets refer to the numbered references listed at the end of 


this Documentation File. 
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D2.2 Zenith Angle at Solar Noon, sn 


The zenith angle at solar noon, ésy, is the angle at solar noon between 
the beam radiation and the normal to a horizontal surface. At solar noon, the 
zenith angle is a function of only the declination, 6, and latitude, ¢. 


Cos Oc, = Cos & Cos ¢ + Sin & Sin ¢ (D-2) 
{see Ref. 24, pg. 13, Eq (1.6.4)] 


D2.3 Sunset Hour Angle, » 


The sunset hour arqle, w, is the angular displacement of the sun, west of 
the local meridian, at sunset due to the rotation of the earth on its axis at 15° 


per hour. 
Cos w = -Tan ¢ Tan & (D-3) 


{see Ref. 24, pg. 13, Eq (1.6.7)] 


02.4 Atmospheric Transmittance for Beam Solar Insolation at Solar Noon, 7p SN 


The atmospheric transmittance for beam solar insolation at solar noon, 
7B,SN» may be determined from Eq (D-10). The constants ag’, aj*, and k” must 
first be determined from Eqs (D-4) through (D-6) for the specified altitude, A. 
The constants rg, rj and ry may be selected from Table D1. 


ay. = 0.4237 - 0.00821 [6.0 - (a/3281)]* / (D4) 
{see Ref. 24, pg. 62, Eq (2.8.2)] 
a," = 0.5055 + 0.00595 [6.5 - (A/3281)]° (D-5) 


{see Ref. 24, pg. 62, Eq (2.8.3)] 


th) 
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k= 0.2711 + 0.01858 [2.5 - (A/3281)]* (D-6) 


{see Ref. 24, pg. 62, Eq (2.8.4)] 


Table D1 - Correction Factors for Climate Type” 


Climate Type 


Tropical 


Mid-Latitude Summer 
Subarctic Summer 


_Mid-Latitude Winter 


* see Ref. 24, pg. 62, Table 2.8.1. 


ap. = ag ro [see Ref. 24, pg. 63] (D-7) 
aj = ay ry [see Ref. 24, pg. 63] | (D-8) 
k= k*ry [see Ref. 24, pg. 63] (D-9) 
7B,SN = a0 + aj exp (-k/Cos &sn) (D-10) 


{see Ref. 24, pg. 62, Eq (2.8.1)] 


02.5 Atmospheric Transmittance for Diffuse Solar Insolation at 
Solar Noon, 7p SN 


The atmospheric transmittance for diffuse solar insolation, at solar noon, 
TD,SN» may be determined by Eq (D-11). 


7D,SN = 0.2710 - 0.29397D, SN (D-11) 


[see Ref. 24, pg. 64, Eq (2.8.7)] 


Copyright por American Petroleum Institute 
Wed Feb 09 12:31:01 2005 


APT MPMS*19.3D) 93 MM 0732290 0511533 30T 


tio of Hourly to Dail Total Solar Insolation on a Horizontal Surface 


D2.6 R y y 


at Solar Noon, ry SN 


The ratio ry,sn is the ratio of the hourly total solar insolation on a 
horizontal surface at solar noon, Iy,sn, to the daily total solar insolation on a 
horizontal surface, Hy. This ratio may be calculated from Eq (0-14), where the 
coefficients a and b are first determined from Eqs (D-12) and (0-13), 
respectively. 


a = 0.4090 + 0.5016 Sin (w - 60) (0-12) 
[see Ref. 24, pg. 79, Eq (2.13.2a)] 
b = 0.6609 - 0.4767 Sin (w - 60) (D-13) 


[see Ref. 24, pg. 79, Eq (2.13.2b)] 


tf (a + b)(1-Cos w) 
TH sn * 5 —_—_————— (D-14) 
i 24| |Sin w - (2xw/360) Cos w 


{see Ref. 24, pg. 79, Eq (2.13.1)] 


02.7 Hourly Total Solar Insolation on a Horizontal Surface. at Solar Noon, IH SN 


The hourly total solar insolation on a horizontal surface at solar noon, 
Iy,SN» may be determined from the ratio TH,SN and the daily total solar 
insolation on a horizontal surface, Hy, using Eq (D-15). 


TH,SN* HSN HH | (D-15) 


02.8 Diffuse Solar Insolation at Solar Noon, In sn 


The diffuse solar insolation at solar noon, Ip,sn» may be determined from 
Iq SN» 7B,SN and 7D, SN using Eq (D-16). 
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: 
DSN 
t a 2.) -= 
Tp,sn ™ 4H,sN se , | oe) 
73,5N ~ "DSN 


D2.9 Beam Solar Insolation at Solar Noon, Ig SN 


The beam solar insolation at solar noon, Ig,SN» may be determined from 
In,sn> Ip,sn and @sn using Eq (D-17). 


I - I 

I ~ |HLSN DSN (0-17) 

8,SN : 
Cos Bon 


D3.0 SAMPLE CALCULATIONS 


This section presents sample calculations illustrating how Eqs (D-1) through 
(D-17) may be used to determine the solar insolation parameters for a sample 
problem. 


The Given Conditions chosen for this sample problem are the same as those 
that were used to develop the API Computer Data Base [39]. 


Given Conditions: 


n = 105 year day number 

= 35 deg. north latitude 
A = 2000 ft. altitude 
Hy = 1500 B/day ft? 


The Calculated Results of this sample problem are summarized in Table D2. 
These results are also used in Section C6 to develop an equation for the vapor 
space temperature factor, Ky, for the case of a “two color fixed-roof tank with 


typical solar insolation parameters" (see Eqs (C-84) through (C-88)). 
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03.1 Calculate § from Eq (0-1) 


360 (284 + 105) 
6 = (23.45) Sin |——— 
365 


= 9.41490° 
D3.2 Calculate @sy from Eq (D-2 


Cos @sy = Cos (9.41490) Cos (35) + Sin (9.41490) Sin (35) 


= 0.90195 
@sy =: 25..5848° 


03.3 Calculate w from Eq (D-3) 
Cos w = - Tan (35) Tan (9.41490) 


= -0.116106 
@ = 96.6674° 


D3.4 Calculate an*, ay” _and k™ from Eqs (D-4) through (D-6) 


ap” = 0.4237 - 0.00821 [6.0 - (2000/3281)]* 


= 0.185144 

aj” = 0.5055 + 0.00595 [6.5 - (2000/3281)]* 
= 0.711948 

K* = 0.2711 + 0.01858 [2-5 - (2000/3281)]2 
= 0.337499 


03.5 Determine rp, ry and rx from Table Dl 


ro = 0.97 
r1 = 0.99 
ry = 1.02 


DS 
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D3.6 Calculate aq, ay and k from Eqs (D-7) through (D-9 


ag = (0.185144) (0.97) 
= 0.179590 

a] = (0.711948) (0.99) 
= 0.704829 

k = (0.337499) (1.02) 
= 0.344249 


03.7 Calculate 7 from Eq (D-10 


7B, SN = 0.179590 + 0.704829 exp (-0.344249/Cos (25.5848) ) 
= 0.660790 


D3.8 Calculate rp sn from Eq (D-11) 


7D,SN = (0.2710) - (0.2939) (0.660790) 
= 0.0767938 


D3.9 Calculate a and b from Eqs (D-12) and (0-13) 


a = 0.4090 + 0.5016 Sin (96.6674 - 60) 
= 0.708540 

b = 0.6609 - 0.4767 Sin (96.6674 - 60) 
= 0.376230 


03.10 Calculate ry sn from Eq (D-14) 


- [(0.708540) + (0.376230)][1 - Cos (96.6674)] 
"H, SN F Sin (96.6674) - (2x)(96.6674/360) Cos (96.6674) 


= 0.133277 
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D3.11 Calculate Ip sn from Eq (0-15) 


TH,SN = (0.133277) (1500) 
= 199.915 B/hr ft2 


D3.12 Calculate Ip sn from Eq (0-16) 


0.0767938 
= (199.915) |———-_..-.... 
0.660790 + 0.07679 


Ip SH 


= 20.8142 B/hr ft? 


D3.13 Calculate Ip syn from Eq (D-17) 


(199.915 - 20.8142) 
B, SN Cos (25.5848) 


= 198.571 B/hr ft? 


D4.0 CONCLUSION 


The typical solar insolation parameters developed in this Section D of the 
Documentation File were used in Section C6 to develop the equation for 
calculating the daily vapor temperature range, ATy. 
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Table 02 - Summary of Sample Problem Results 


Symbol Description 


Given Information 


Year Day Number 
Latitude deg. norgh 
Daily Total Solar Insolation on a Horizontal Surface B/day ft 
Altitude ft 


Calculated Results 


Declination Angle — 9.41490 
Zenith Angle at Solar Noon 25.5848 
Sunset Angle : 


oa 0.344249 


Atmospheric Transmittance for Beam Solar Insolation 
at Solar Noon 0.660790 
Atmospheric Transmittance for Diffuse Solar Insolation 


at Solar Noon 0.0767938 
mae 0.708540 


--- 0.376230 
Ratio of Ip sy to Hy 0.133277 
Hourly Total Solar Insolation on a Horizontal Surface 

at Solar Noon 199.915 B/hr ft2 
Diffuse Solar Insolation at Solar Noon 20.8142 B/hr ft2 
Beam Solar Insolation at Solar Noon 198.571 B/hr ft2 
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Fp Paint factor dimensionless 
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p Solar reflectance of paint dimensionless 
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E1.0 INTRODUCTION 


This section of the Documentation File to API Publication 2518, Second 
Edition, contains a development of the solar absorptance, a, for selected paint 
colors, paint types and paint conditions. The values developed are based on 
information that was presented in Appendix IV of API Publication 2518, First 


Edition [A6]*. 


Section E2 presents fundamental definitions and a relationship between solar 
absorptance and solar reflectance. 


Section €3 describes the work that resulted in the paint factors, Fp» that 
appear in API Publication 2518, First Edition. 


Section £4 describes the relationship between solar absorptance and paint 


factor. 


Section £5 develops a relationship between the solar absorptance of paints 
in poor condition and good condition. 


Section EE summarizes the set of solar absorptance values that were selected 
for use in API Publication 2518, Second Edition. 


£2.0 SOLAR ABSORPTANCE FUNDAMENTALS 


Solar radiation that impinges on a surface is either: (1) absorbed by the 
surface; (2) reflected from the surface; or (3) transmitted through the surface. 
The fraction of the incident solar radiation (referred to as "insolation”) that 
is absorbed is referred to as the absorptivity, a, of the surface; the fraction 
that is reflected is referred to as the reflectivity, p, of the surface; and 
the fraction that is transmitted is referred to as the transmittance, 7, of the 
surface, where: 


Leen aan 


. Numbers in brackets refer to the numbered references listed at the end of 


this Documentation File. 
£4 
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atptre=l (E-1) 
For the metal surfaces of fixed-roof tanks, the transmittance is zero, or: 


r=0 (£-2) 
Eq (E-1) then reduces to: 
atp=l (E-3) 


The following discussion focuses primarily on the values of absorptance, 
a, of a surface for use in calculating the thermal breathing losses of fixed- 
roof tanks. If the reflectance, p, of a surface is known, Eq (E-3) may be used 
to determine the corresponding absorptance, a, of the surface. 


The exterior surface of fixed-roof tanks is normally coated with a paint to 
reduce corrosion. A wide range of paint colors have been used, sometimes with a 
different color on the tank roof than on the tank shell. 


The absorptance of tank paint depends upon the paint color, paint type and 
paint condition. Newly painted tank surfaces, or painted surfaces in a good 
condition will have a lower absorptance than weathered painted surfaces or 
painted surfaces in poor condition. 


E3.0 DEVELOPMENT OF PAINT FACTORS FOR API PUBLICATION 2518, FIRST EDITION 


At the time that the First Edition of API Publication 2518 {A6] was 
published, the importance of the effect of paint absorptance on the thermal 
breathing loss was recognized. A paint with a low absorptivity, such as white 
paint, was known to affect the thermal breathing loss in two significant ways: 


1. It reduces the transfer of heat to and from the tank vapor space and 
therefore reduces the volume of thermal breathing loss. 


2. It reduces the transfer of heat to the bulk liquid and therefore 
reduces the thermal breathing loss by lowering the stock vapor 
pressure. 
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During the development of the First Edition to API Publication 2518, a Paint 
Factor Task Group was formed and was given the assignment of developing a 
practical correlation of the effect of paint on the thermal breathing loss from 
fixed-roof tanks. Through discussions and correspondence, the task group worked 
with paint chemists and the staff of one large paint manufacturer. In addition, 
pertinent literature was studied and analyzed. Evaporation loss tests, other 
than those used in developing the breathing loss correlation [see Eq (3) in Ref. 
A6], were studied and tabulated [see Table IV-2 in Ref. AG]. Tests were 
performed to determine the difference in liquid body temperature in large tanks 
painted aluminum versus those painted white. Tests were performed with both 
artificial light and sunlight on various painted surfaces to determine their 


reflectance. 


In this discussion, reflectance refers to measurements made with a- 
spectrophotometer which utilizes an integrating sphere. On this scale, 100 
represents the reflectance of pure magnesium oxide (MgO) sprayed on tile until no 
increase in the instrument reading can be observed. This instrument measures 
only the visible spectrum. Although heat for evaporation is partially supplied 
to a tank by the nonvisible (infared) portion of the sun’s rays, it was assumed 
that the effect of this portion is proportional to the visible (measured) portion 
of the sun’s rays.. It was not concluded, however, that this relationship will 
hold for all paint formulations and climatic conditions. 


In developing the correlation for thermal breathing loss [Eq (3) in Ref. 
A6], test data from 64 painted tanks were selected for detailed analysis. The 


paint classification for the 64 tanks tested is as follows: 


Table El - Paint Classification for Tanks Tested 


Number 
of Tanks 


White White 
Aluminum White 
White Aluminum 
Aluminum Aluminum 
White Gray 
Gray Gray 


Total 
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At the time that the First Edition of API Publication 2518 was published, little 
loss data on fixed-roof tanks with paint colors other than white, aluminum or 
gray was available. 


Early in the study of paint reflectance, a definite relationship appeared to 
exist between relative loss and tank paint reflectance. It also became obvious 
that the broad paint classifications of white, aluminum, gray, etc. were not 
sufficiently definitive. Paints have a wide range of reflectance, varying from a 
freshly painted laboratory sample to a surface in a dirty and extremely weathered 
condition. 


Paints of the same nominal color and in the same condition have a wide range 
of reflectance, depending upon formulation. One paint manufacturer reported that 
the reflectance of new white paint on a tank might range between 0.80 and 0.86, 
with a maximum of 0.88, and that the reflectance of new black paint on a tank 
might range between 0.04 and 0.08. New aluminum paint also has a wide range of 
reflectance, depending upon the grade of aluminum powder and the type of 
formulation. One type of aluminum paint, formulated with a specular: type of 
pigment (polished scales), may have a reflectance that might range between 0.60 
to 0.70, with an average of 0.68. Another type of aluminum paint, utilizing a 
diffuse type of pigment, may have a reflectance as low as 0.35. Aluminum paints 
of intermediate reflectance can be formulated from mixtures of these two paint 
types. Thus, aluminum color paints cover a wide range on the total reflectance 
scale. 


Approximate reflectance values for various paint colors and paint conditions 
are presented in Table E2. These values are based on: (1) the thermal breathing 
loss correlation test data [see Tables II-1, [1-2 and I1I-3 in Ref. A6]; (2) 
supplementary test data; (3) information supplied by paint manufacturers; and, 
(4) the results of field tests conducted by the task group. 


For reference purposes, Table E3 lists the reflectance values of various 
paint colors and paint conditions reported by Nelson {3]. 


Table £4 summarizes the set of paint factors, Fp, that were developed by the 
Paint Factor Task Group for use with the thermal breathing loss correlation [Eq 


£7 


Copyright por American Petroleum Institute 
Wed Feb 09 12:31:04 2005 


API MPMS*19.1D) 93 @@ 0732290 O5S1154b Tbd 


(3) in Ref. A6). These paint factors were compiled by the task group from a 
judicious review of all available data. These paint factors do not represent a 
precise evaluation of the effect of tank paints on evaporation loss from tanks. 
They do, however, present useful and reasonable paint factors that are based on 
the alignment of the relative test data and on the agreement of the test data for 
aluminum painted tanks and gray-painted tanks. 


£4.0 RELATIONSHIP BETWEEN SOLAR ABSORPTANCE AND PAINT FACTORS 


Figure €1 {from Fig. IV-3 of Ref. AG}. illustrates the effect of solar 
reflectance, p, on the paint factor, Fp. This figure shows a_ linear 
relationship between p and Fp for values of Fp < 1.5. Eq (E-3) expresses 
this linear relationship: 


Fp = 1.74 - 0.90 p (E-3) 


Substituting Eq (E-2) into (E£-3) and solving for the solar absorptance, a, we 
obtain: 


@ = 1.11 Fp - 0.94 (E-4) 


Eq (E-4) was used to convert the values of paint factor, Fp, to solar 
absorptance, a, in Table £4. 


£5.0 EFFECT OF PAINT CONDITION ON SOLAR ABSORPTANCE 


The solar absorptance values listed in Table £4 include values for paint in 
good condition, ag, and values for paint in poor condition, ap. 


Figure £2 is a plot of ap versus ag for each of the paints listed in 
Table &4. The values plotted may be fit with a linear relationship that must go 
through the point where ag = 1 and ep = 1. Eq (E-5) expresses this linear 
relationship: 


ap = 0.2 + 0.8 ag (E-5) 
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E6.0 CONCLUSION 


Table ES summarizes the set of solar absorptance values that were selected 
for use in API Publication 2518, Second Edition. The values selected were based 
upon a careful evaluation of the values presented in Tables £2 and £4. 


Use was made of the values in Table E4 for only the cases where the tank 
roof and shell are painted the same color. Also, Eq (E-5) was used to determine 
the poor condition values from the good condition values for gray and red color 


paints. 


The values of Nelson {3] listed in Table £3 were not incorporated into Table 
ES, but are listed in this section of the Documentation File only for reference. 
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Table E2 - Solar Reflectance and Solar Absorptance of Tank Paints (A6] 


Solar 
Absorptance 


Solar 
Reflectance 


Condition 
p (1) @ (2) 
(dimensionless) (dimensionless) 


1.00(1) 


Magnesium 
Oxide (3) 


White 


0.83 
0.65 


Aluminum Specular Good 


0.60 
0.45 


Specular Average 


Aluminum 


Diffuse Good 


Aluminum 


0.40 


Diffuse Average 


Aluminum 


Extremely 0.30 


Weathered 


Diffuse 


Aluminum 


0.11 
0.06 


Red Primer Good 


Black 


Notes: (1) Values of p are from Table IV-1 of Ref. A6. 
(2) Values of @ are calculated from the relationship: 
a = l-p 
(3) Magnesium oxide is the reflectance standard whose value is 
assumed to be 1.00. 
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Table E3 - Solar Reflectance and Solar Absorptance of Tanks Paints {3] 


Selar 
Reflectance 


Solar 
Absorptance 


Paint 
Condition 


a p 
(dimensionless) (dimensionless) 


Aluminum New 


Aluminum 
Aluminum 


Black 


Blue 
Blue 
Blue 
Blue 


Weathered 


Light 
Pale 
Dark 


Cream Light 


Glossy 
Light 


Gray 
Gray 
Gray 


Light 


Green 
Green 
Green 
Green 


Dark 
Dark 


Bare 


Metal 


Pink Light 


Red Iron Oxide 
Red 
Red 


Red 


Dark 
Bright 


Tan 
White 


Yellow 
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Table £5 - Solar Absorptance for Selected Tank Paint Colors and Conditions 


Solar Absorptance, a (dimensionless) 
Paint Condition 


Aluminum | Specular 


Aluminum | Diffuse 


Light 


Gray 


Gray Medium 


Red Primer 


White 
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i 

Ze 
Li 
tiiA7At 


T | Vit 
Het 


100 0.80 0,60 0,40 0.20 0 


Point Factor, Fp (dimensionless) 


Solar Reflectance, » (dimensionless) 


Paint Factor 
i El - Effect of Solar Reflectance on 
eae [From Fig. IV-3 of Ref. A6] 


E14 
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Solar Absorptance of Poor Condition Paint, 
ap (dimensionless) 


Solar Absorptance of Good Condition Paint, 
ag (dimensionless) 


Figure £2 - Effect of Paint Condition on Solar Absorptance 
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SECTION F 
NOMENCLATURE 
SYMBOL DESCRIPTION UNITS 
A Area ft2 
a Defined by Eq (F-21) dimensionless 
b Defined by &q (F-22) — dimensionless 
¢c Heat capacity B/1bm °F 
G Defined by Eq (F-9) dimensionless 
Hy Daily total solar insolation on a horizontal 
surface B/day ft2 
h Heat transfer coefficient B/hr ft2 OF 
KL Defined by Eq (F-6) dimensionless 
k Thermal conductivity B/hr ft °F 
T Temperature OF 
AT Temperature change OF 
Solar absorptance dimensionless 
p Density lbm/ft3 
SUBSCRIPTS 
A Ambient 
AA Ambient, Daily Average 
B Liquid Bulk 
BA Liquid Bulk, Daily Average 
BN Liquid Bulk, Daily Minimum 
BX Liquid Bulk, Daily Maximum 
G Gas 
GA Gas, Daily Average 
GN Gas, Daily Minimum 
GX Gas, Daily Maximum 
L Liquid Surface 
LA Liquid Surface, Daily Average 
LN Liquid Surface, Daily Minimum 
LX Liquid Surface, Daily Maximum 
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F1.0 INTRODUCTION 


This section of the Documentation File to API Publication 2518, Second 
Edition, presents the development of the equations required to determine the 
daily maximum liquid surface temperature, TLx, the daily average liquid surface 
temperature, TLa, and the daily minimum liquid surface temperature, TLN- 


Section F2 develops an equation for estimating the liquid bulk temperature, 
Tg, when this temperature is not available from tank operating records. This 
| equation is based on data from in API Publication 2518, First Edition (A6]*. 


Section F3 presents an equation for determining the daily average gas space 
temperature, TéA- 


Section F4 develops an equation for determining the daily liquid surface 
temperature range, AT_- This theoretical equation is confirmed by 4 
correlation developed from the API Computer Data Base [39]. 


Section F5 develops an equation for determining the daily average liquid 
surface temperature, Tia. This theoretical equations is confirmed “by the API 
Test Data [33,38]. 


Section F6 presents a summary of the equations required to calculate the 
daily maximus liquid surface temperature, Try, the daily average liquid surface 
temperature, TLA, and the daily minimum liquid surface temperature, TiN: The 
equations are presented in the order required for calculation. 


F2.0 LIQUID BULK TEMPERATURE 


-The liquid bulk temperature, Tg, is the average temperature of the liquid 
stock in the storage tank. This information is usually available from tank 
gaging records or other tank operating records. 


TS 


. Numbers in brackets refer to the numbered references listed at the end of 


this Documentation File. 
F4 
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If the liquid bulk temperature is not available, it may be estimated from 
the daily average ambient temperature, Tag, and the tank shell paint solar 
absorptance, a, using Eq (F-1). 


Tg = Ta + 62 - 1 | (F-1) 


Figure Fl is a plot of liquid bulk temperature, Tg, above that in a white 
tank versus shell paint solar absorptance, a. This figure is from Figure IV-2 
of API Publication 2518, First Edition [A6]. Eq (F-1) was developed from the 
straight line on Figure Fl and from the assumption that the liquid bulk 
temperature in a white tank is the same as the average ambient temperature, TAA- 


F3.0 AVERAGE GAS SPACE TEMPERATURE 
The following equation may be used to determine the daily average gas space 


temperature, Tea. Refer to &q (C-108) in Section C8 of this Documentation File 
for the development of this equation. 


Tog = 0.775T ag + 0.225TLA + 0.0140 @ Hy. 


F4.0 LIQUID SURFACE TEMPERATURE RANGE 


(F-2) 


This section contains the development of two equations for the daily liquid 
surface temperature range, ATL$ (1) a correlation from the API Computer Data 
Base [39]; and (2) a theoretical equation from a heat transfer analysis. 


F4.1 Liquid Surface Temperature Range Correlation 


A correlation between the liquid surface temperature range, ATL, and the 
gas space temperature range, ATg, was developed from the 560 sets of calculated 
results that are contained in the API Computer Data Base [39]. 


The daily liquid surface temperature range, ATL, and the daily gas space 
temperature range, ATg, are defined by Eqs (F-3) and (F-4). 
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AT, = TLX - TUN (F-3) 


ATg = Tox - TEN (F-4) 

Figures F2, F3 and F4 are plots of liquid surface temperature range, ATL , 
versus gas space temperature range, ATg, from the API Computer Data Base [39]. 
Fach of these figures contains about 1/3 of the 560 calculated values in the API 
Computer Data Base. 


The resulting correlation is shown in Eq (F-5): 


AT, = 0.5340504Tg - 2.810841 (F-5) 


Although the data in Figures F2 through F4 exhibit a slight curvature, the linear 
correlation results in an excellent fit of the values, with a correlation 


coefficient of 0.998344. 


F4.2 Liquid Surface Temperature Range Equation 


This section contains the development of a theoretical equation for the 
liquid surface temperature range, ATL, from a heat transfer analysis. 


The API Test Data [33,38] indicates that the top layer of liquid stock 
behaves thermally as if it were a layer of solid material during the daily 
thermal heating cycle. Heat is transferred in a vertical direction through this 
top layer by steady periodic heat transfer. 


The case of heat transfer between a fluid medium with a steady periodic 
temperature change in contact with the plane surface of a semi-infinite solid 
material is treated by Jakob [Jakob, M., "Heat Transfer", Volume 1, John Wiley & 
Sons, Inc., 1949, pp. 296-299]. This analysis may be applied to the top layer of 
liquid stock in a storage tank and results in &q (F-6), (F-7) and (F-8): 
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F-6) 
Ky = mas . ( 
AT, 
or 
F-7 
ATL = Ky Alg (F-7) 
where: 
1/2 . 
K : 1 ; (F-8) 
tL lhe 26/7? 4 26 
¢ F at “ (F-9) 
= ey ae 
24 he 


The following values for pL, CL and k_ from the API Base Case [39] were used to 
calculated G and K,: 


pL = 53 Iba/ft? (F-10) 
cL = 0.45 Btu/lbm °F (F-11) 
ky = 0.08 Btu/hr ft OF (F-12) 
hg = 0.65 Btu/hr ft2 OF | (F-13) 


The resulting values calculated for G and K, from Eqs (F-9) and (F-8), 
respectively, are: 


F-14 
G = 0.5911 (F-14) 


-15 
Ky = 0.5185 (F-15) 
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This value for K, is supported by the correlation result of Eq (F-5). 


The value for Ki listed in Eq (F-15) may be rounded to 0.500 for ease of 
calculation. Finally, substituting this value for Ki into Eq (F-7), we obtain: 


AT = 0.500 ATg (F-16) 


F5.0 AVERAGE LIQUID SURFACE TEMPERATURE 


FS5.1 Heat Transfer Model 


This section develops a simplified equation for estimating the daily average 
liquid surface temperature, Tia, from the liquid bulk temperature, Tg, and the 
daily average gas space temperature, Tga- 


Figure F5 is a schematic of the energy flows and temperatures near the 
liquid surface. The heat transfer equations are Eqs (F-17) and (F-18). 


ag = he Ay (Te - Th) (F-17) 
au = hy Ay (Ty - TB) (F-18) 

The heat balance equation at the liquid surface is: 
(F-19) 


aq = 4 
Substituting Eqs (F-17) and (F-18) into Eq (F-19) and solving for Ty we obtain: 


TL = a Tg + b Tg (F-20) 
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he 
bz | | (F-22) 
ha +h 


ace (F-23) 


Substituting Eq (F-23) into Eq "(F-22) we obtain the following relationship 
between b and a. 


b=l-a (F-24) 


During the daily thermal breathing process, the liquid bulk temperature, Tz, 
can be assumed to be constant. At the maximum, average and minimum gas space 
temperature conditions, Eq (F-20) becomes: 


TLx = aTgx + bTgB (F-25) 
TLA = alga + bTB (F-26) 
TLN = algy + bTg (F-27) 


The daily liquid surface temperature range, ATL, and the daily gas space 
temperature range, ATg, were defined by Eqs (F-3) and (F-4). Substituting Eqs 
(F-25) and (F-27) into Eq (F-3), we obtain: 

AT, = (aTgx + bTg) - (aTgn + bTp) 


AT, = a(Tex - Ten) 


AT, = adTg (F-28) 
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From Eq (F-16) in Section F4.2, we know that a = 0.500. Substituting this 
value for a into Eq (F-28), we obtain b = 0.500. Substituting these values for a 
and b into Eqs (F-25), (F-26) and (F-27), we obtain: 


Tix = 0.500Tex + 0.500Tg (F-29) 
TLA = 0.500Tg, + 0.500TR (F-30) 
TLN = 0.500Tgy + 0.500Tg (F-31) 


F5.2 Comparison With Test Data 


_ To compare the results of this analysis with the API Test Data [33,38], it 
is first convenient to form the ratio (TLa - Tg)/(TGa - Tp) using Eq (F-30) from 
the above analysis as follows: 


27 (0.500T,, + 0.500T,) - T 
fa _,| eran #928001 p I Te) 5 k59 _(F-32) 
Tea > T Toa = Tp 


Table Fl presents the above ratio calculated from API Tests 1 through 6. 
API Tests 7 through 10 could not be used fcr this comparison since the liquid 
surface thermocouple appeared to be reading the gas temperature rather than the 
liquid surface temperature. The average value calculated for the ratio (TLa - 
Tg)/(Tea - Tp) from the API Test Data is 0.467. This value calculated from test 
data is close to the value of 0.500 in Eq (F-32), thus verifying the heat 
transfer analysis in Section F5.1. 


F5.3 Average Liquid Surface Temperature Equation 


Substituting Eq (F-2) for Tega into €q (F-30) and solving for Tia, we obtain 
finally the following expression for the daily average liquid surface 
temperature, Tia. 


TLA = 0.437Tag + 0.563Tg + 0.00789aH} (F-33) 
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F6.0 CONCLUSION 


The equations that are required to ca 


below. These equations are presented in the order that is 


calculation. 
1. Determine Tp 


2. 


- Tg = value from tank operating data 


or 
Tg = Tan + 6a - 1 


Determine TLA 


TLA = 0.437Tpa + 0.563Tg + 0.00789aH, 


Determine ATy 


ATy = 0.723AT, + 0.0279aH} 


Determine ATj 
ATy = 0.500ATy 
Determine Ty 


Tix = TLA + 0-500ATy 


Determine Tin 


Tin = TLA - 0-500ATL 


Tculate Tra, TLx and TLN are summarized 


required for 


(F-34) 


(F-35) 


(F-36) 


(F-37) 


(F-38) 


(F-39) 


It is possible to combine some of the above equations to reduce the number 
of equations, but they are listed separately here for clarity. 
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Table Fl - Temperature Difference Ratio (TLA - TBA)/(TGA - TBA) 


Toa-TBA 


(dimensionless) 


Calculated From the API Test Data [33,38] 


API-1 
API-2 
API-3 
API-4 
API-5 
API-6 


Average 


Note: (1) The mean values Tga, TLA» Tya listed in this table were 
calculated from the maximum and minimum values measured 


in the API tests [33,38]. 
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(Liquid Bulk Temperature - Average Ambient Temperature), (Tg-Taa)(°F) 


0 0.2 0.4 0.6 0.8 1.0 


Solar Absorptance, a (dimensionless) 


Figure Fl - Effect of Paint Solar Absorptance on the 
Liquid Bulk Temperature [A6] 
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Figure F5 - Schematic of Energy Flows and 
Temperatures Near the Liquid Surface 
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SENSITIVITY ANALYSIS OF 
STANDING STORAGE LOSS EQUATION 


Gl 
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SECTION G 
NOMENCLATURE 
SYMBOL DESCRIPTION 
_A Constant in the vapor pressure equation 
B Constant in the vapor pressure equation 
D Tank diameter 
Ey Daily standing storage loss calculated from 
API Publication 2518, First Edition [A6] 
E9 Daily standing storage loss calculated from 
API Publication 2518, Second Edition [A7] 
Hyo Vapor space outage 
I Daily total solar insolation on a horizontal 
surface 
Ke Vapor space volume expansion factor 
Ks Vented vapor saturation factor 
Ls Daily standing storage loss 
My Stock vapor molecular weight 
PA Atmospheric pressure 
Pep Breather vent pressure setting (always a 
positive value) 
Pay Breather vent vacuum setting (always a 
negative value) 
Pya Stock vapor pressure at the daily average 
liquid surface temperature 
Pyn Stock vapor pressure at the daily minimum 
liquid surface temperature 
Pyy Stock vapor pressure at the daily maximum 
liquid surface temperature 
APy Stock daily vapor pressure range 
APg Breather vent pressure setting range 
RVP Reid vapor pressure 
TAA Daily average ambient temperature 
TB Liquid bulk temperature 
TLA Daily average liquid surface temperature 
ATA Daily ambient temperature range 
TiN Daily minimum liquid surface temperature 
G3 
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UNITS 


dimensionless 
OR 
ft 


lb/day 


1b/day 
ft 


Btu/ft2 day 
dimensionless 
dimensionless 
]b/day 
1b/1bmole 
psia 


psig 
psig 
psia 
psia 


psia 

psi 

psi 

psi 

OR or OF 
OR or OF 
OR or OF 
OR or OF 
OR or OF 


API MPMS*19.1D 93 MM 0732290 0531574 922 mw 


API PUBLICATION 2518 
DOCUMENTATION FILE 


| SECTION G 
| 
| NOMENCLATURE (Continued) 
SYMBOL DESCRIPTION 
TLX Daily maximum liquid surface temperature 
Aly Daily vapor temperature range 
Vy Tank vapor space volume 
Wy Stock vapor density 
a Tank paint solar absorptance 


UNITS 


OR or OF 

OR or OF 

ft3 

1b/Ft3 
dimensionless 
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G1.0 INTRODUCTION 


This section of the Documentation File to API Publication 2518, Second 
Edition [A7]}* contains a sensitivity analysis of the standing storage loss 
equation. The objective of this analysis is to determine the affect of the 
variables in the standing storage loss equation on the calculated loss over 4 
range of conditions for each variable. 


Section G2 summarizes the variables used in the sensitivity analysis. 


— 


Section 63 presents sample calculations for the Base Case conditons. 


SS 


Section G4 presents the results of the sensitivity analysis. 


———— TT 


Section G5 presents conclusions from the sensitivity analysis. 


62.0 VARIABLES USED IN THE SENSITIVITY ANALYSIS 


Table Gi lists the input variables that were studied as part of the 
sensitivity analysis. For each input variable, a "Base Case" value and a range 
were selected. 


In performing the sensitivity analysis, each of the input variables was kept 
at the Base Case value while the value of the input variable being studied was 
varied over the desired range. In this way, the sensitivity of the API standing 
storage loss model was evaluated around the Base Case conditions. 


63.0 SAMPLE CALCULATIONS FOR THE BASE CASE — 


This section summarizes sample calculations of the standing storage loss for 
the Base Case conditions. 


oe ee ee 
{ = Numbers in brackets refer to the numbered references listed at the end of 


this Documentation File. 
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63.1 Base Case Conditions 


63.1.1 Tank Conditions 


D = 100 ft 
Hyg = 30 ft 
a = 0.50 


Pgp = 0.03 psig 
Peay = -0.03 psig 


63.1.2 Meteorological Conditions 
Tag = 519.69°R (60.00°F) 
ATa = 20.00°R (20.00°F) 
I = 1500 Btu/ft2 day 
Pa = 14.696 psia 
63.1.3 Stock Conditions 
Type = Crude Oi] 
RVP = 2.0 psi 
My = 80 1b/lbmole 


63.2 Calculated Results 


63.2.1 Determine Tg 


Tg = Tag + 6a - 1 
Tg = (519.69) + (6)(0.50) - (1) 
Tg = 521.69°R (62.00°F) 
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63.2.2 Determine TiA 
TLA = 0.44Tag + 0.56Tg + 0.0079al 
Tia = (0.44)(519.69) + (0.56) (521.69) + (0.0079) (0.50) (1500) 
Tia = 526.74°R (67.05°F) 
63.2.3 Determine Aly 
ATy = 0.72ATa + 0.028al 
aTy = (0.72) (20.00) + (0.028) (0.50) (1500) 
ATy = 35.40°R (35.40°F) 
63.2.4 Determine TLx 
TLX = TLA + 0.25ATy 
Try = (526.74) + (0.25) (35.40) 
TLy = 535.599R (75.90°F) 
63.2.5 Determine TiN 
Tin = TLA - 0.25ATy 
Tin = (526.74) - (0.25)(35.40) 
TiN = 517.89°R (58.20°F) 
63.2.6 Determine A 
A = 12.82 - 0.9672 In (RVP) 


A = (12.82) - (0.9672) In(2.0) 
A = 12.15 


G7 
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63.2.7 Determine B 


B = 7261 - 1216 In(RVP) (G-7) 
B = (7261) - (1216) In(2.0) 
B = 6418°R 
63.2.8 Determine Pyx 
: G-8) 
Pyy = exp [A - * ( 
LX | 
(6418) 
P,. = exp|(12.15) - ———— 
vx = oP) (535.59) 
= ].1817 psia 
63.2.9 Determine Pya 
B 
G-9 
Pua = exP A - =a (6-9) 
LA 
(6418) 
Py, = exp! (12.15) - ———— 
va = xy (526.74) 


= 0.9662 psia 
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63.2.10 Determine Pyn 


A (6-9) 
= exp|A - — 


(6418) 
= exp} (12.15) - ——— 


: ~ (517.89) 


VN 


= 0.7846 psia 


63.2.11 Determine APy 


(6-11) 
APy = Pyx - PyN 


APy = (1.1817) - (0.7846) 
APy = 0.3971 psi 


G3.2.12 Determine APp 


(G-12) 
APg = Pep - Pay 

APp = (0.03) - (-0.03) 

APg = 0.06 psi 


G3.2.13 Determine Vy 


c (G-13) 
2 


Wg 
(3.1416) (100)2(30) 
(4) 


= 235,620 ft? 


Vy = 


Vy 


G9 
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t 


G3.2.14 -Determine Wy 


My Pua (G-14) 
RTp 
(80.00) (0.9662) 


iy = ———_—_— 
(10.731) (526.74) 
~ 0.01367 1b/Ft? 


V 


Wy 


63.2.15 Determine Ke 


Tra Pa” Pv 
(35.40) (0.3971) - (0.0600) 
E (526.74) (14.696) - (0.9662) 
Kr = 0.09176 
63.2.16 Determine Ks 
Ke = ee eee (G-16) 
1 + 0.053 Py, Hyg 
l 
1 + (0.053) (0.9662) (30) 
K, = 0.3943 
63.2.17 Determine Ls 
Ls = Vy Wy Ke Ks (G-17) 


ls = (235, 620) (0.01367) (0.09176) (0.3943) 
Ls = 116.5 1b/day 
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64.0 SENSITIVITY ANALYSIS RESULTS 


Table G2 lists the input and output variables from the API standing storage 
loss equation for each of the 57 cases which were studied as part of the 
sensitivity analysis. 


Figures Gl through G10 illustrate the effect of the input variables on the 
emissions calculated by the API standing storage loss equation. The Base Case 
condition is also indicated on each figure. 


64.1 Breather Vent Setting Sensitivity 


Figure 68 indicates that the breather vent pressure and vacuum settings (up 
to t0.05 psig) have little effect on the standing storage loss. 


G4.2 Vapor Pressure Sensitivity 


The sensitivity analysis also include an evaluation of the effect of vapor 
pressure on emissions. Vapor pressure is affected by stock type, Reid Vapor 
Pressure and temperature. Table G2 and Figures G9 and G10 show the effect of 
Reid Vapor Pressure and resulting true vapor pressure, respectively, on emissions 
for crude oil. 


Even though the RVP was kept constant in most of the sensitivity analysis, 
the TVP varied from 0.3547 psia in Case No. 29 to 2.2815 psia in Case No. 33. 
This variation in TVP occurred implicitly as the calculated liquid surface 
temperature vaired due to varying the ambient temperature. 


Currently, stocks with a TVP greater than 1.5 psia are not commonly stored 
in fixed-roof tanks. 
64.3 moarison Between the First and Second Editions of API Publication 2518 


Comp 


A comparison was made between standing storage loss calculated from the 
First Edition [A6] and the Second Edition [A7] of API Publication 2518. The RVP 


G11 
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was varied from 0.001 psi to 12.0 psi, while the other variables were kept at the 


Base Case conditions. 


Table G3 summarizes the calculated results. The last column indicates the 
ratio of the standing storage loss calculated from the First Edition to that 
calculated from the Second Edition. It should be noted that the ratio equals 
about 1.0 when the RVP equals about 12 psi or the TVP equals about 10 psia. 


Figures G9 and G10 graphically illustrate the effect of RVP and TVP on 
standing storage loss. It should be noted that the emissions predicted by the 
Second Edition are close to the emissions predicted by the First Edition when the 
vapor pressure is over about 2 psia. For low vapor pressures, the emissions 
calculated from the Second Edition are significantly less than those calculated 
from the First Edition, indicating that the First Edition overpredicts emissions 
for low vapor pressure stocks (less than 2.0 psia). 


65.0 CONCLUSIONS 


The sensitivity analysis showed that the standing storage loss equation 
exhibits the same trends that were shown by a similar sensitivity analysis 
performed on the API Computer Model (39]. This similarily in sensitivity helps: 
validate the standing storage loss equation. 


The sensitivity analysis showed that the breather vent pressure and vacuum 
settings had little affect on the calculated loss. 


The sensitivity analysis showed that for high vapor pressure stocks (greater 
than 2.0 psia) the evaporation losses calculated from the Second Edition are 
close to those calculated from the First Edition, and that for low vapor pressure 
stocks (less than 2.0 psia) the First Edition overpredictes the evaporation 


loss. 
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Table 63 - Effect of Reid Vapor Pressure on Emissions Calculated From 
From Both the First and Second Edition of API Publication 2518 


Ej Ee 
(5) (6) 
(1b/day) (1b/day) 

15,660 | 0.00002705 | 0.00003607 0.00759 
14,320 | 0.0001223 | 0.0001591 0.03312 
12,860 | 0.0006232 | 0.0007893 0.1664 
11,520 | 0.002817 | 0.003481 0.7249 
10,060 | 0.01450 0.01745 3.599 
8,725 | 0.06428 0.07546 14.89 
7,261 | 0.3335 | 0.3811 59.67 
6,418 | 0.8587 0.9662 116.5 
5,925 | 1.496 1.668 168.9 
5,304 | 2.983 3.289 282.1 
4,732 | 5.694 6.212 534.5 
4,461 | 7.682 8.339 846.4 


. « « ry ry e . ry ry 
S o o So So to) o w fend 


4,239 | 9.919 10.723 1529. 


Notes: (1) Calculated from: A = 12.82 - 0.9672 In (RVP). 
(2) Calculated from: B = 7261 - 1216 In (RVP). 
(3) Calculated from: Pg = exp [A - (8/TB)» where Tg = 521.69°R (62.00°F). 
(4) Calculated from: Pya = exp [A - (B/TLq)], where TLa = 526.74 (67 .05°F). 
(5) Calculated from API Publication 2518, First Edition [A6} for the 
Base Case Conditions listed in Table Gl and for the RVP listed in 
this table. 
(6) Calculated from API Publication 2518, Second Edition [A7] for the 
err re Conditions listed in Table Gl and for the RVP listed in 
this table. 
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Tank Diameter, D (ft) 


Figure Gl - Effect of Tank Diameter, J, on Emissions 
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Figure G2 - Effect of Vapor Space Outage, Hyp, on Emissions 
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Figure G5 - Effect of Stock Vapor Molecular Weight, My, on Emissions 
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API PUBLICATION 2518 
DOCUMENTATION FILE 


SECTION H 


COMPARISON OF STANDING STORAGE LOSS 
EQUATION WITH TEST DATA 


NOTE: This document is part of the API standards development process and is 
intended for use by API committees ONLY. It shall not be reproduced or 
circulated or quoted, in whole or in part, outside of the cognizant API 
committee(s) except with the written approval of API. This draft API standard 
will be formatted and edited prior to API publication. Copyright @1990. 
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SECTION H 
NOMENCLATURE 
SYMBOL DESCRIPTION 
A Constant in the vapor pressure equation 
B Constant in the vapor pressure equation 
0 Tank diameter 
Hyo Vapor space outage 
I Daily total solar insolation on a horizontal 
surface . 
Ke Vapor space expansion factor 
Ls Daily standing storage loss 
My Stock vapor molecular weight 
Pa Atmospheric pressure 
Pgp Breather vent pressure setting (always a 
positive value) 
Pay Breather vent vacuum setting (always a 
negative value) 
Q Vented gas volume outflow 
Tax Daily maximum ambient temperature 
TAN Daily minimum ambient temperature 
Tg Stock liquid bulk temperature 
Aly Daily vapor temperature range 
Vy Tank vapor space volume 
a Tank point solar absorptance 
H3 


Copyright por American Petroleum Institute 
Wed Feb 09 12:31:17 2005 


UNITS 
dimensionless 
OR 
ft 
ft 


Btu/ft2 day 
dimensionless 
lb/day 
1b/1bmole 
psia 


psig 


psig 

ft3/day 

OR 

OR 

OR 

OR 

ft3 
dimensionless 
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H1.0 INTRODUCTION 


This section of the Documentation File to API Publication 2518, Second 
Edition [A7]* contains a comparison of the standing storage loss equation with 
the API test data [38], EPA test data [20] and WOGA test data [17]. 


H2.0 COMPARISON WITH THE API TEST DATA [38] 


Table Hl summarizes the input information from Tests API-1 through API-10 
which was used to calculate the standing storage loss. This table also compares 
the measured and calculated results. 


All of the 10 API tests [38] continued sufficient information to make a 
comparison between the measured and calculated results for ATy, Q and Ls. . 


Table Hl summarizes a comparison between the measured and calculated vapor 
space temperature range, Aly. The average difference is -17.5%. Figure Hl 
illustrates this comparison graphically. 


Table Hl summarizes a comparison between the measured and calculated vented 
gas volume outifow, Q, which was calculated from Q = Vy Ke. The average 
difference is -11.3%. Figure H2 illustrates this comparison graphically. 


Table Hl summarizes a comparison between the measured and calculated daily 
standing storage loss, Ls. The average difference is -14.3%. Figure 4H3 


illustrates this comparison graphically. 


In general, the comparison between the API test data and the predictions 
made by the standing storage loss equation are excellent. 


a 


= Numbers in brackets refer to the numbered references listed at the end of 


this Documentation File. 
H4 
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H3.0 COMPARISON WITH THE EPA TEST DATA {20} 


Table H2 summarizes the input information from the EPA tests which was used 
to calculate the standing storage loss. This table also compares the measured 
and calculated results. 


Only 12 of the 15 EPA tests [20] were suitable for making a comparison 
between the measured and calculated results for ATy, Q and Ls. 


Table Hl summarizes a comparison between the measured and calculated vapor 
space temperature range, ATy. The average difference is 10.7%. Figure H4 
sTlustrates this comparison graphically. 


Table Hl summarizes a comparison between the measured and calculated vented 
gas volume outflow, Q. The average difference is 77.3%. Figure H5 illustrates 
this comparison graphically. 


Table Hl summarizes a comparison between the measured and calculated daily 
standing storage loss, Ls. The average difference is 111.3%. Figure 4H6 
illustrates this comparison graphically. 


In general, the thermal response of the tank vapor space is predicted quite 
well, thus indicating that the heat transfer aspects of the standing storage loss 
model compare well with actual measurements. 


Hé.0 COMPARISON WITH THE WOGA TEST DATA [17] 


Table H3 summarizes the input information from the WOGA tests which was used 
to calculate the standing storage loss. This table also compares the measured 
and calculated results. 


During the WOGA tests, measurements were not made of the vapor space 


temperature, Ty. Thus, it is not possible to compare the predicted values with 
the measured values for this variable. 


H5 
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Only 8 of the 44 WOGA tests {17] were suitable for making a comparison 
between the measured and calculated results for Q and Ls. 


Table H3 summarizes a comparison between the measured and calculated vented 
gas volume outflow, Q. The average difference is 79.3%. Figure H7 illustrates 


this comparison graphically. 


Table H3 summarizes a comparison between the measured and calculated daily 
standing storage loss, Ls. The average difference is 28.1%. Figure H8 
illustrates this comparison graphically. 


In general, the comparison between the measured and calculated vented gas 
volume outflow, Q, for the WOGA tests is about as good as that for the EPA tests. 
Again, this illustrates that the standing- storage loss model incorporates a good 
charaterization of the tank vapor space thermal response. 


H5.0 CONCLUSIONS 


The comparison with the API test data [38], EPA test data [20] and WOGA 
test data [17] validated the suitability of the standing storage loss equation. 


The API tests provided a more accurate and extensive set of test data for 
comparison with the API standing storage loss equation. The average difference 
between the calculated and measured standing storage loss is 14.3%. 


The EPA and WOGA test data, although of lesser suitability for an accurate 
validation, also confirmed the suitability of the standing storage loss equation. 
The average difference between the calculated and measured standing storage loss 
is. 111.3% for the EPA tests and 28.1 for the WOGA tests. 
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Calculated Versus Measured Vapor Space Temperature Range, 
for the API Tests [38] 


ATy, 


10 “20 3 
Measured Vapor Space Temperature Range, ATy, (OF) 


Figure Hl - 


Oo i=) cS oS 
wo = on N one! 


(Jo) ‘Ayy ‘abuey aunjzeusdwa, aseds sodea peze_nayed 


H10 


Copyright por American Petroleum Institute 
Wed Feb 09 12:31:19 2005 


API MPMS¥19.1) 93 0732290 0511607 1T7 = 


oa 
wo 
NI 


Measured Vented Gas Volume Outflow, Q (ft3/day) 
for the API Tests [38] 
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Measured Daily Standing Storage Loss, Ls (1b/day) 


Figure H3 - Calculated Versus Measured Daily Standing Storage Loss, Lc, 
For The API Tests [38] 
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Figure H4 - Calculated Versus Measured Vapor Space Temperature Range, Ty, 


for the EPA Tests [20] 


APT MPNS419.1) 93 MM 0732290 05116410 7953 


Measured Vented Gas Volume Outflow, Q (ft3/day) 


Figure HS - Calculated Versus Measured Vented Gas Volume Outflow, Q, 
for the EPA Tests [20] . 
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Measured Daily Standing Storage Loss, Ls (1b/day) 
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Measured Daily Standing Storage Loss, Ls .(1b/day) 


Figure H8 - Calculated Versus Measured Daily Standing Storage Loss, Ls, 
for the WOGA Tests [17] 
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R1 INTRODUCTION 


This section of the Documentation File to API: Publication 2518, Second 
Edition, contains a list of the references that are cited in this Documentation 
File. In addition to those that are cited, the list also contains references 
seis were reviewed as part of preparing the Second Edition to API Publication 
2518. 


Section R2 contains a list of API publications that relate to atmospheric 
storage tanks. These are arranged in numerical order according to their API 
publication number. 

Section R3 contains a list of literature references. These are arranged in 
chronological order according to their publication date. 


R2 AMERICAN PETROLEUM INSTITUTE PUBLICATIONS 


Al. American Petroleum Institute, "Inspection_of Atmos heric and _Low-Pressure 
Storage Tanks", Recommended Practice RP-575, First Edition, Washington, 
D.C., 1990. ° 


AZ. American Petroleum Institute, “Welded Steel Tanks for Oi] Storage", Standard 
650, Seventh Edition, Washington, D.C., November 1980. 


A3. American Petroleum Institute, "Venting Atmospheric and Low-Pressure Storage 
Tanks (Non refrigerated and Refrigerated)", Standard 2000, Third Edition, 
Washington, 0.C., January 1982, Reaffirmed December 1987. 

Aé. American Petroleum Institute, "Evaporation Loss in the Petroleum Industry- 
Causes and Control", Bulletin 2513, First Edition, Washington, 0.C., 

ebruary 1959. 

AS. American Petroleum Institute, "Evaporative Loss from External Floating-Roof 
Tanks", Publication 2517, Third Edition, Washington, D.C., February 1989. 

AS. American Petroleum Institute, “Eva oration Loss from Fixed-Roof. Tanks", 
dis 2518, First Edition, Washington, D.C., June 1962, Reaffirmed August 

A7. American Petroleum Institute, "Evaporative Loss from Fixed-Roof Tanks”, 
Publication 2518, Second Edition, Washington, 0.c., vecember 1990. 


AS. American Petroleum Institute, “Evaporation Loss from Internal Fioating-Roof 
Tanks", Publication 2519, Third Edition, Washington, D.C., June 1983. 


AQ. American Petroleum Institute, “Use of Pressure-Vacuum Vent Valves for 
Atmospheric Pressure Tanks to Reduce Evaporation Loss", Bulletin 2521, First 


Edition, Washington, D.C., September 1966. 


Al0. American Petroleum Institute, "Petrochemical Evaporation Loss from Storage 
Tanks", Bulletin 2523, First Edition, Washington, D.C., November 1969. 
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All. American Petroleum Institute, "Guide for Inspection of Refinery Equipment", 
Chapter XIII, Atmospheric and Low-Pressure Storage Tanks, Fourth Edition, 


10. 


11. 


Washington, D.C., April 1981. 

LITERATURE REFERENCES 

Boardman, H.C., “Storage of Volatile Petroleum Products", Petroleum Refiner, 
Vol. 25, No. 4, April 1946, pp. 109-116. 

American Petroleum Institute, Symposium on Evaporation Loss of Petroleum 
from Storage Tanks, Papers Preseated During the 32nd Annual Meeting of the 


American Petroleum Institute, Held in Chicago, IL, November 10, 1952, 
(Also Published in API Proceedings, Vol. 32, Part I, 1952, pp. 212-261). 


Nelson, W.L., “How Painting Affects Storage Tank Losses”, The Oi1 and .Gas 
Journal, Vol. 52, November 2, 1953, pg. 130. 


Prater, N.H., “How to Calculate Vapor Losses", Petroleum Processing, April 


1954, pp. 537-540. 


Walker, E.H., R.M. Eltringham, and A. Puttick, "Evaporation Loss From Petrol 
Storage Tanks in the United Kingdom - A Practical Survey", Paper Presented 
at a Meeting of The Institute of Petroleum at 26 Portiand Place, London, 
U.K., May 11, 1955. 

Hoffman, E.L., “The Effect of Paint Colors in Reducing Storage Tank Losses”, 
Paper Presented at the 35th Annual Meeting of the American petroleum 
Institute, Marketing Division, Operations and Engineering Committee, San 
Francisco, California, November 14, 1955. 


Snyder, A.D., “Theoretical Analysis of Breathing Volatile Liquid Vapor 
Losses from Fixed-Roof storage Tanks", Report, Prepared by Olin Mathison 
Chemicals Co. for the American Petroleum Institute, Committee on Evaporation 
Loss, July 7, 1965. 


Cinnamon, S.J., “Investigation of a Method for Reducing Vapor Losses _in 
Storage of Volatile Fiquids”. N.S. Thesis, Submitted to Purdue University, 
Department of Chemical Engineering, August 1965. 

Cinnamon, S.J. and J.E. Myers, “Reducin Liquid Storage Losses", Chemical 
Engineering Progress, Vol. 61, No. 12, Decenber 1965, pp. 69-74. 

Air Pollution Control Association, “Control of Atmos heric Emissions from 
Petroleum Storage Tanks", Informative Report No. 2, Prepared by the TI-3 
Petroleum Committee, Presented in the Journal of the Air Pollution Control 
Association, Vol. 21, No. 5, May 1971, pp. 260-268. 

NUKEM GmbH, “Development of a Method for the Recognition of Environmental 
Damage Using the Example of Emissions Due to Transfer of Fuels and Breathin 

of Tank Tretallations”. Prepared for the Ministry of the Interior, Federal 


a 
Republic of Germany, 1972. 


Copyright por American Petroleum Institute 
Wed Feb 09 12:31:22 2005 


12. 


13. 


14. 


15. 


16. 


17. 


18. 


19. 


20. 


21. 


22. 


23. 


API MPMS*¥19-1) 93 MM 0732290 0511618 T34 


Matsumura, I., “Evaporation Loss of Hydrocarbons in _Handlin Petroleum", 
Bulletin of the Japan Petroleum Institute, Vol. 16, No. 2, November 1974, 


pp. 132-139. 


BMI-DGMK, “Measurement and Determination of Hydrocarbon Emissions in the 
Course of Storage and lransfer in Above-Ground Fixed Cover Tanks With and 
Without Floating Covers, Project 4590-11: Fixed Cover tank Without Floatin 

Cover®, jointly funded by BMI (Ministry o the Interior, Federal Republic o 

Germany) and DGMK (German Society for Petroleum Technology and Coal 
Chemistry), 1976. 


Ball, D.A., A.A. Putman and R.6. Luce, “Evaluation of Methods for Measuring 
and Controlling Hydrocarbon Emissions from Petroleum Storage Janks", 
Battelle Columbus Laboratories, Prepared for the U.S. Environmental 


Protection Agency, EPA Report No. 450/3-76-036, November 1976. 


Zanker, A., “Estimate Tank Breathing Loss", Hydrocarbon Processing, January — 
1977, pp. 117-120, (Errata contained in Hydrocarbon Processing, July 1977, 
pp. 82-83). 


Schwanecke, R., “Evaporation Losses from Liquid Tanks", Wasser, Luft and 
Betrieb, Vol. 15, No. 7, July 1977, pp. 254-258. 

Engineering-Science, Inc., "H drocarbon Emissions From Fixed-Roof Petroleum 
Tanks", Prepared for the Western Oil and Gas Association, July 1977. 

Harrer, R.D., "Field Data Developed for Fixed-Roof Tank Emissions", The Oil 
and Gas Journal, January 2, 1978, pp. 90-97. 

Wilson, A.L., "Suggested Emission Factors for Fixed-Roof Storage Tanks", 
Paper Presented at the Meeting of the Air Pollution Control Association, 


Anaheim, California, November 13, 1978. 

Engineering-Science, Inc., "Synthetic Organic Chemical Manufacturing 
Industry, Emission Test Report, Breathin Loss Emissions from Fixed-Roof 
Petrochemical Storage Tanks", Prepared for the U.S. nvironmental Protection 
Agency, EPA Report No. MB-78-0CM-5, February 1979. 

Moryzkov, V.S., L.L. Tatarnikov, E.Y. Kardash and A.S. Yarmukhametov, 
"tosses of Crude Oi] and Products in Operation of Refinery Storage Tanks", 
Khimiya 1 Tekhnologiya Topliv i Masel, No. 4, April 1979, pp. 8-10. 


Cinquemani, V., J.R. Qwenby, Jr., and R.G. Baldwin, “Input for Solar 
Systems”, Prepared by the U.S. Dept. of Commerce, National Oceanic and 
Atmospheric Administration, Environmental and Information Service, National 
Climatic Center, Asheville, North Carolina, Prepared for the U.S. Dept. of 
Energy, Div. of Solar Technology, under Interagency Agreement No. E (49-26)- 
1041, November 1978 (Revised August 1979). 


Erbar, J.H., “Predicting Hydrocarbon Emissions from Fixed Roof Field Storage 
Tanks", Prepared for the American Petroleum Institute, Environmental Affairs 
Department, Task Force on Hydrocarbon Emissions from Production Operations, 
1980. 
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Duffie, J.A. and W.A. Beckman, "Solar Engineering of Thermal Processes”, 
John Wiley and Sons, 1980. 


Laverman, R.J., “Evaluation of the Temperature Factor for Thermal Breathing 
Loss Estimation from Fixed Roof Tanks’, icago Bridge & Iron Co., CBI 
Research Report No. RC-3000-2, March 11, 1981. 

Laverman, R.J., "Thermal Breathing Loss Analysis for a Roof Vent Area 
Diffusion Rate Controlled Model", Chicago Bridge & Iron Co., CBI Research 
Report No. RC3000-1, March II, 1981. 

Beckman, J.R. and J.R. Gilmer, "Model for Predicting Emissions From Fixed- 
Roof Storage Tanks", Industrial & Engineering Chemistry, Process Design & 
Development, Vol. 20, No. 4, April 1981, pp. 646-651. : 

U.S. Environmental Protection Agency, “Compilation of Air Pollutant Emission 
Factors, Section 4.3, Storage of Organic Liquids”, Supplement No. 12 to EPA 
eport AP-42, April 1981. 

TRW Environmental, Inc., “Background Documentation for Storage of Organic 
Liquids", Prepared for the U.S. Environmental Protection Agency, EPA 
Contract No. 68-02-3174, May 1981. 

Beckman, Duffie and Associates, “Evaporation Loss of Petroleum from Storage 
Tanks", Final Report, Prepared for the American Petroleum Institute, 


Committee on Evaporation Loss Measurement, August 1, 1982. 


ites Deutscher Ingenieure, “Emission Control: Refineries", VDI 2440, June 


Beckman, J.R., "Breathing Losses from Fixed-Roof Tanks by Heat and Mass 
Transfer Diffusion, Industrial & Engineering Chemistry, Process Design & 
Development, Vol. 23, No. 3, 1984, pp. 472-479. 


Environmental Monitoring & Services, Inc. (subsidiary of Combustion 
Engineering Co.), "Breathing Loss Emissions from Fixed Roof Tanks", Final 
Report, Prepared for the American Petroleum Institute, Committee on 
Evaporation Loss Measurement, June 1985. 

U.S. Environmental Protection Agency, “Compilation of Air Pollutant Emission 
Factors, Section 4.3, Storage of Organic Liquids”, USEPA Report No. AP-42, 
Third Edition, September 1985. 

U.S. Dept. of Commerce, National Oceanic and Atmospheric Administration, 
"Comparative Climatic Data Through 1984", National Climatic Data Center, 
Asheville, North Carolina, 1986. 


Beckman, J.R. and J.A. Holcomb, “Experimental and Theoretical Investigation 
of Working Emissions from Fixed-Roof tanks", industrial and Engineering 
Chemistry, Process Design & Development, Vol. 25, No. 1, 1986, pp. 293-298. 


Copyright por American Petroleum Institute 
Wed Feb 09 12:31:22 2005 


37. 


38. 


39. 


40. 


APT MPMS*1L9-1D 93 MM 0732290 0511620 b32 


Beckman, J.R., "Model Development to Predict Hydrocarbon Emissions from 

Crude Oi] Storage and Treatment Tanks", Final Report, Prepared for the 

Ea onus Air Resources Board, CARB Agreement No. A4-045-32, July 20, 
1986. 

Knodel, B.D. and R.J. Laverman, “Data Base Generation, Analysis and Revision 
f API Publication 2518, Task 1: Validate Computer Model*, Final Report for 


) 
Task 1, Prepared by CBI Industries, Inc., Prepared for the American 


- Petroleum Institute, Committee on Evaporation Loss Measurement, Task Group 


2518, September 11, 1986. 


Rinehart, J.K. and R.J. Laverman, "Data Base Generation, Analysis and 
Revision of API Publication 2518, Task 2: Generate Computer Data Base”, 
Final Report for task 2, Prepared by CBI Industries, Inc., Prepared for the 
American Petroleum Institute, Committee on Evaporation Loss Measurement, 
Task Group 2518, February 16, 1987. 


Rinehart, J.K. and R.J. Laverman, “Data Base Generation, Analysis and 
Revision of API Publication 2518, Task 3: Correlate Data Base", Final Report 

or task 3, Prepared by CBI Industries, Inc., Prepared for the American 
Petroleum Institute, Committee on Evaporation Loss Measurement, Task Group 


2518, August 26, 1988. 
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